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We consider a minimal scale-invariant extension of the Standard Model of particle physics com- 
bined with Unimodular Gravity formulated in [1]. This theory is able to describe not only an 
T— H inflationary stage, related to the Standard Model Higgs field, but also a late period of Dark Energy 

^~H domination, associated with an almost massless dilaton. A number of parameters can be fixed by 

^^ inflationary physics, allowing to make specific predictions for any subsequent period. In particular, 

Cn we derive a relation between the tilt of the primordial spectrum of scalar fluctuations, n^, and the 

1—^ present value of the equation of state parameter of dark energy, lJt,^. We flnd bounds for the scalar 

^ ? tilt, ns < 0.97, the associated running, —0.0006 < dlnua/dlnk < —0.00015, and for the scalar-to- 

' tensor ratio, 0.0009 < r < 0.0033, which will be critically tested by the results of the Planck mission. 

T-H For the equation of state of dark energy, the model predicts lJoe > ~1- The relation between Us 

y~~^ and ojqb allows us to use the current observational bounds on n^ to further constrain the dark 

energy equation of state to < 1 -I- a;^^; < 0.02, which is to be confronted with future dark energy 
' _H ' surveys. 
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I. INTRODUCTION 



At the classical level, the Lagrangian describing the Standard Model of particle physics (SM) minimally 
coupled to General Relativity (GR) contains three dimensional parameters: Newton's constant G, the vacuum 
expectation value (vev) of the Higgs field or, equivalently, the Higgs boson mass and a possible cosmological 
[^ constant A. The masses of quarks, leptons and intermediate vector bosons are induced by the vev of the 

Higgs field. At the quantum level, additional scales, such as Aqcd and all other scales related to the running 
of coupling constants, appear due to dimensional transmutation. It is tempting to look for models in which 
all these seemingly unrelated scales have a common origin. 

In this work we present a detailed analysis of a model realizing this idea, proposed in [1]. We will refer 
to it as the Higgs-Dilaton model. The model is based on a minimal extension of the SM and GR that 
; ^ contains no dimensional parameters in the action and is therefore scale-invariant at the classical level. Scale 

Cd invariance is achieved by introducing a new scalar degree of freedom, called Dilaton. The motivation of the 

model relies on the assumption that the structure of the theory is not changed at the quantum level. In 
other words, the full quantum effective action should still be scale-invariant and the effective scalar potential 
should preserve the features of the classical potential. A perturbative quantization procedure maintaining 
scale invariance was presented in [2] (see also [3]). In the Higgs-Dilaton model all scales are induced by 
the spontaneous breakdown of scale invariance (SI). As a consequence of the broken symmetry, the physical 
dilaton is exactly massless. Replacing GR by Unimodular Gravity (UG), in which the metric determinant 
is fixed to one, \g\ — 1, results in the appearence of an arbitrary integration constant in the equations of 
motion, representing an additional breaking of scale symmetry. As discussed in [1], in theories with scalar 
fields non-minimally coupled to gravity, this constant effectively gives rise to a non-trivial potential for the 
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scalar fields. In the case of the Higgs-Dilaton model, the new potential is of the "run-away" type in the 
direction of the dilaton. 

While the dynamical breakdown of the scale symmetry by the Higgs field can provide a mechanism for 
inflation in the early universe [4], the light dilaton, practically decoupled from all SM fields, can act as 
Quintessence (QE), i.e. as dynamical Dark Energy (DE). We find that, under some assumptions, it is 
possible to relate the observables associated to infiation to those associated to dark energy. Namely, we 
establish a functional relation between the predicted value for the tilt n^ of the primordial scalar power 
spectrum and the predicted equation of state parameter w^^ of dark energy. Further, we find a relation 
involving the corresponding second order quantities, i.e the running a^ of the spectral tilt and the rate of 
change wf)^ of the DE equation of state. 

The paper is organized as follows: In section II A we introduce and discuss the minimal scale-invariant 
extension of the Standard Model and General Relativity. In Section II B the idea of Unimodular Gravity is 
described and applied to the scale-invariant model. We then discuss the cosmology of the resulting Higgs- 
Dilaton Model. The infiationary period is studied in detail in Section III. The implications of the model for 
the late dark energy dominated stage are studied in Section IV. Finally, conclusions are presented in Section 
V. For completness, an analysis of slow-roll infiation in the Jordan frame and the differences with respect 
to the Einstein frame are presented in Appendix A. 

II. THE HIGGS-DILATON MODEL 

In this section we review the Higgs-Dilaton model of [1], which consists of two moderate extensions of the 
Standard Model and General Relativity (SM plus GR) . In subsection II A we show how the introduction of a 
dilaton allows to extend SM plus GR to a phenomenologically viable scale-invariant theory. After discussing 
the main properties of the resulting theory (II A 1), we discuss two naturalness issues, the Cosmological Con- 
stant Problem and the Gauge Hierarchy Problem, in the context of this model (II A 2). Next, we give some 
arguments in favor of a particular parameter choice corresponding to the absence of a cosmological constant 
(II A 3). In subsection II B the construction of the model is completed by replacing GR by Unimodular 
Gravity (UG). The qualitative picture of cosmology in the Higgs-Dilaton model, as found in [1], is recalled 
in subsection II C. 



A. Minimal scale-invariant extension of SM plus GR 

1. Introducing the Dilaton 

Let us start by writing down the Lagrangian density that combines GR and the SM-'^ 

-^ = \mIR + £sM[A^o] - A (^V -«')'- A , (1) 

where the first term is the usual Einstein-Hilbert action for GR with Mp = (SttG)"'^/^, the second term is 
the SM Lagrangian without the Higgs potential, the third term is the Higgs potential with the SM Higgs 
doublet Lp and its vacuum expectation value v, and A is a cosmological constant. In this standard theory, to 
which we will refer as "SM plus GR", classical scale invariance is violated by the presence of the dimensional 
constants M , v and A. Our goal is to let these scales be dynamical, i.e. replace them by a field. The most 
obvious solution, without introducing new degrees of freedom, would be to let the Higgs field be responsible 
for all scales. This corresponds to considering the Lagrangian^ 
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£,tp^ipR + LsM[\^o] - A ((^V) , (2) 
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1 We use the conventions r;^^ = diag{-l, 1, 1, 1) and _R"„ ^ = d^Vi^ + V^ V\g - (7 ^ &). 
^ With the conventions used here, a conformally coupled scalar field has ^ = —1/6. 



£_ being a new real parameter ("non-minimal coupling"). The associated action is now scale-invariant, i.e. 
invariant under the global transformations, 

gp^(x) H^ g^uicrx) , $(x) H^ a'^'^^iax) , (3) 

where <I>(x) stands for the different particle physics fields, d^ is their associated scaling dimension and a 
is an arbitrary real parameter. In a theory that is invariant under all diffeomorphisms (Diff invariant), as 
is the case for (2), the symmetry associated to the absence of dimensional parameters can equivalently be 
written as an internal transformation^ 

gt^uix) ^ a-^g^,^{x) , $(x) ^ cr'^*$(x) . (4) 

Can the Lagrangian (2) give a satisfactory phenomenology? Since we are looking for a theory that should 
eventually be quantized, we want to introduce the requirement that the theory possesses a "Classical Ground 
State". The term "Classical Ground State" will be used throughout this work to refer to solutions of the 
classical equations of motion, which correspond to constant fields in the particle physics sector of the theory 
and a maximally symmetric geometry, i.e. Minkowski (flat), de Sitter (dS) or Anti de Sitter (AdS) spacetime. 
The existence of such a ground state might be essential for a consistent quantization of the theory. At the 
quantum level, the theory should possess a ground state that breaks scale invariance and in this way induces 
masses and dimensional couplings for the excitations (particles). We will require that this spontaneous 
symmetry breaking already appears in the classical theory due to the existence of a symmetry breaking 
classical ground state'*. 

Let us now look for symmetry-breaking classical ground states in the theory (2). If gravity is neglected, 
i.e. the first term in the Lagrangian is dropped, the classical ground states correspond to the minima of 
the scalar potential A (f^ip) ■ The only possibility for them to break the scale symmetry, ip = ipo ^ 0, 
is to set A = 0. In this case the theory possesses an infinite family of classical ground states satisfying 
2(p^ip = Hq, where h^ is an arbitrary real constant. If one includes gravity, the set of possible classical 
ground states becomes richer. Namely, even if A 7^ the theory possesses a continuous family of classical 
ground states satisfying 2(p'''ip = Hq and R = 4A/iq/^, where ho is an arbitrary real constant. The states 
with ho ^ break scale invariance spontaneously and induce all scales at the classical level. Hence, the 
goal of having a classical theory in which all scales have the same origin, spontaneous breakdown of SI, 
is achieved. However, the above theory is in confiict with experimental constraints. In fact, although the 
non-zero background value of (p gives masses to all other SM particles, the excitations of the Higgs field itself 
are massless and, moreover, decoupled from the SM fields [7]. This fact is seen most easily if the Lagrangian 
is written in the Einstein- frame by defining the new metric ^^^ = Mp (.f^ifigfj^i, and the new canonical Higgs 
field (p = Mp\Jl/£, + 6 ln((^ /Mp). (This type of variable change will be discussed in detail in section III.) 
In the new variables, the SI of the original formulation corresponds to a shift symmetry for the Higgs field 
(p, which is the massless Goldstone boson associated to the spontaneous breakdown of SI. A Higgs field with 
these properties is excluded by Electroweak precision tests [8]. Therefore, in order to construct a viable SI 
theory, it seems unavoidable to introduce new degrees of freedom. 

The next simplest possibility is to add a new singlet scalar field to the theory. We will refer to it as the 
dilaton x- The scale-invariant extension for the SM plus GR including the dilaton reads 

4= ^\ (^xx' + 2a^V) R + /:sM[A^o] - ^g^^a^xa.x - v{p, x) , (5) 

where the scalar potential is given by^ 

n^,x) = A(^V-|^x')'+/3x^ (6) 



^ Note that in a theory that is invariant only under a restricted class of diffeomorphisms, such as UG, the absence of dimensional 

parameters will still guarantee invariance under (3) but not under (4). 
^ The authors of [5, 6] propose that scale symmetry could be broken by the pure presence of a time-dependent cosmological 

background. 
^ The parametrization chosen for the scalar potential assumes A 7^ 0. This only excludes the phenomenologically unacceptable 

case where a quartic term (ip^ip)^ is absent. 



We will only consider positive values for ^^ and ^h, such that the coefficient in front of the scalar curvature is 
positive, whatever values the scalar fields take. The positivity of the non-minimal coupling parameters is at 
the same time the condition for semi-positive definitness of the scalar field kinetic terms. By construction, 
the action associated to (5) is invariant under (3) , respectively (4) . The theory should possess a symmetry- 
breaking classical ground state with ip ^ ipo ^ and X — Xo 7^ 0. The case ipo = would correspond to 
a theory with no electroweak symmetry breaking, while the case Xo = would result in a theory with a 
massless Higgs field. Both these cases are phenomenologically unacceptable. 

Let us again start by neglecting the gravitational part of the action. In its absence, the ground states 
correspond to the minima of the potential (6). It is easy to see that the only possibility to get a ground 
state satisfying (po ^ and xo 7^ is to have a potential with a flat direction, i.e. a > and f3 — 0, as well 
as A > for stability. The corresponding family of classical ground states is given by 2(f^ip = h^ and x = Xo 
with Hq = ^Xo' where xo is an arbitrary real constant. 

Like before, the inclusion of gravity results in the appearance of additional classical ground states for 
/3 7^ 0, given by 

A A A4x + aS,h 

The solutions with Xo 7^ spontaneously break SI. All scales are induced and proportional to xo- For 
instance, one can directly identify the Planck scale as 

f2 _ i ,2 , <- , 2 _ /^<- , . " , 4/3eI \ ,2 



Depending on the value of /3, the background corresponds to flat spacetime (/3 — 0), de Sitter or Anti de 
Sitter spacetime of constant scalar curvature R, corresponding to a cosmological constant 






A ^ Uil,R ^ ^^^^ _ . (9) 



The spectrum of perturbations around a symmetry-breaking solution contains the usual massless spin-2 
perturbation in the gravitational sector. The scalar sector contains an excitation with mass 

which will play the role of the physical SM Higgs field, plus a massless Goldstone boson (both perturbations 
are combinations of the fields x and h). We use h to denote the field ip in the unitary gauge. Like in 
the standard Higgs mechanism, the excitations of the Standard Model fields get masses proportional to Hq. 
If one extends the SM by introducing right-handed neutrinos [9, 10], these neutrinos get induced masses 
proportional to xo^- 

In the described model, physics is completely independent of the value of xoi as long as xo 7^ 0. This 
is because only dimensionless ratios of the different scales can be measured. Therefore, parameters of the 
model have to be chosen such that these ratios correspond to the measured ones. For instance, one should 
reproduce the hierarchies between the cosmological scale and the electroweak scale, i.e. A/m* ^ 0(10^'"'®), 
as well as the ratio between the electroweak and the gravitational scale m? /Mp ~ 0(10"^^). We choose 
the parameter /3 to be responsible for the first ratio and a for the second ratio. Therefore, these parameters 
have to take values satisfying /3 <Sg: a <gc 1 and /3,a <^ ^x^^h- One then gets approximately A/m'^ ~ 4- 

and 171^ /Mp ~ |^. Note that the order of magnitude relation vA/m^ ^ rn^/Mp respectively -y//3 ^ a is 
the big number coincidence pointed out by Dirac [11]. However, the present model does not address the 
question about the origin of the big differences between theses scales, i.e. the smallness of a and (3, nor does it 
explain their approximate relation. The non-minimal couplings ^^ and ^h will be constrained by cosmological 
considerations, and A < 0(1), as it corresponds to the self coupling of the Higgs field. Therefore, one can 



Gauge invariance does not allow for couplings of x to SM fields. 



fix the values of a and /3 that give the correct ratios. In the same fashion, one has to choose values for the 
SM Yukawa couplings that produce the observed mass ratios. 

As the theory contains a new massless degree of freedom, the dilaton, one has to make sure that it does 
not contradict any experimental bounds. A detailed analysis of the interactions between this massless field 
and the SM fields is contained in [12]. Let us cite the relevant findings of that work. It turns out that as a 
consequence of SI, the massless scalar field completely decouples from all SM fields except for the Higgs field. 
Since the massless field is the Goldstone boson associated to the broken scale symmetry, there exists a set of 
field variables in terms of which it couples to the physical Higgs field only derivatively. In addition, for an 
appropriate choice of field variables, these interactions appear as non-renormalizable operators, suppressed 
by the scale Mp/^h- The analysis of section III will show that S,h ^ 10"''. The suppression scale of non- 
renormalizable operators is therefore considerably lower than the Planck scale, but still much larger than all 
known Particle physics scales. 

Other deviations from the SM appear as a consequence of the non- minimal couplings to gravity. In fact, 
the physical Higgs field, i.e. the field that couples to the SM degrees of freedom, is not h, but a combination 
of h and x- It was shown in [12] that the resulting deviations from the SM are suppressed by the ratio m?' /Mp 
between the physical Higgs mass and the Planck mass, respectively by the small parameter a. While the 
new massless field hardly affects SM phenomenology, we will see that it might play an important role in 
cosmology. 

At the classical level, the above theory successfully implements the idea that all scales are consequences 
of the spontaneous breaking of SI. All conclusions remain true if SI and the features of the potential can 
be maintained at the quantum level (in this context, see ]2, 13, 14]). In that case, the presented model is a 
viable effective field-theory extension of the SM and GR. 

2. Naturalness issues 

The presented theory contains two important fine tunings related to the very big differences between 
the Planck scale Mp, the electroweak scale m and the cosmological scale A. At the quantum level, this 
can lead to two much-discussed naturalness issues. One of them is part of the Cosmological Constant 
Problem. In standard SM plus GR the effective cosmological constant is the sum of a bare constant and 
radiative corrections proportional to the particle physics mass scales of the theory, e.g. the electroweak 
scale. Matching the effective cosmological constant with its observed value, tiny compared to, for instance, 
the electroweak scale, requires a tremendous fine-tuning of the bare cosmological constant. In the case of the 
scale-invariant theory discussed here, the situation is somewhat different. Exact SI forbids a term y/—gA in 
the action. Also, if the quantization procedure respects SI, such a term is not generated radiatively. However, 
as we saw above, due to the non-minimal couplings of the scalar fields to gravity, the cosmological constant 
is in fact associated to the term Px^- Now, this term is not forbidden by scale invariance. Therefore, even if 
scale invariance can be maintained at the quantum level, the quantum effective potential will contain a term 
/3effX^i where /3eff is a combination of the bare value of P and other non-dimensional couplings of the theory. 
These other couplings are generally much bigger than the value of /3eS that corresponds to the observed 
cosmological constant. So, again a strong fine-tuning is needed in order to keep /3eff sufficiently small. This 
tells us that the Cosmological Constant Problem also exists in an exactly scale-invariant theory of the type 
proposed here. 

The second naturalness issue is related to the mass of the Higgs boson and is commonly called "Gauge 
Hierarchy Problem". The problem is twofold. The effective field theory combining the SM with GR contains 
two extremely different mass scales, namely, the electroweak scale v = 246 GcV {v being the vacuum 
expectation value of the Higgs field) and the Planck scale Mp = (STrG)"^/^ = 2.44 • 10^^ GeV. It is 
considered unnatural to have such a huge difference between two scales of the same theory. This is the first 
part of the Gauge Hierarchy Problem. In the considered type of scale-invariant theories, the big difference 
between the electroweak and the Planck scale remains unexplained. 

The other part of the Gauge Hierarchy Problem is related to the stability of the Higgs mass against 
radiative corrections (for a recent discussion see e.g. [15]). Much like the cosmological constant, the mass of 
the Higgs field gets radiative corrections proportional to the other particle physics mass scales of the theory. 
The logic is the same as in the case of the cosmological constant. If there exists a particle physics scale much 
bigger than the electroweak scale, the measured value of the electroweak scale can only be explained by an 
important fine-tuning of parameters. In other words, if there exists a new particle physics scale between 
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the electroweak scale m and the Planck scale Mp, the "smallness" of the Higgs mass constitutes a serious 
theoretical issue. This issue still appears in an exactly scale-invariant theory with spontaneous breaking of 
the scale symmetry. 

If the theory contains no intermediate particle physics scale between m and Mp, the situation is different. 
In that case, whether or not the Higgs mass should be expected to contain big radiative corrections of 
the order Mp depends on the ultraviolet (UV) completion of the theory. At the level of the low-energy 
effective field theory, the UV properties can be encoded in the choice of the renormalization scheme. A 
renormalization scheme based on the assumption that the UV completion is scale-invariant, and which does 
not bring in extra particle physics scales, was presented in [2] (see also [3]). If this scheme is applied to the 
considered minimal scale-invariant extension of SM plus GR, the Higgs mass does not obtain corrections 
proportional to Mp (induced by the vacuum expectation value of the dilaton) and there is no problem of 
stability of the Higgs mass against radiative corrections. Hence, SI makes for the absence of this part of the 
Gauge Hierarchy Problem. 

3. The special case /? = 

We now want to give some arguments in favor of the case jS — 0. This case corresponds to the existence 
of a flat direction in the Jordan- frame potential (6) and hence to the absence of a cosmological constant. 

The reasoning of the precedent paragraph tells us that choosing /3 = corresponds to a fine-tuning of 
the parameters, especially at the quantum level, just like putting A = in standard SM plus GR. From 
this point of view, such a parameter choice should clearly be disfavored. Nevertheless, we think that the 
case /3 = is specially interesting. One reason is that only if /3 = 0, SI can be spontaneously broken in the 
absence of gravity. Put in other words, /3 = allows flat space-time together with [ip, x) = (V'O; Xo) 7^ (0, 0) 
to be a classical solution. 

Another argument is related to the stability of the ground state. As discussed above, a scale-invariant 
theory with spontaneous symmetry breaking always contains a massless scalar degree of freedom, Goldstone 
boson, independently of the value of j3. Now, if /3 7^ 0, the background spacetime of the theory corresponds 
to de Sitter (or anti de Sitter) space-time. It is known, however, that a massless scalar fleld is unstable 
in de Sitter spacetime [16]. There are also indications that this is the case for the 4-dimensional AdS [17]. 
Therefore, it is conceivable that a consistent quantization of the theory might rely on the requirement /3 = 
and hence the existence of flat space-time as a solution (see also ]18-22]). 

A third aspect appears in the context of cosmology. The theory with /3 = 0, not containing a cosmological 
constant, does not seem to withstand the confrontation with cosmological observations. Just like the case 
/? < (AdS) it can not explain the observed accelerated expansion of the universe without introduction of a 
new dark energy component. From this point of view, the only viable option seems to be /3 > (dS). This 
conclusion is correct if gravity is described by GR. However, as we will see in the upcoming section, the 
situation is very different if GR in (5) is replaced by Unimodular Gravity. In that case, the appearance of an 
arbitrary integration constant will give rise to a potential for the Goldstone boson of broken scale invariance. 
As a consequence, for appropriate parameter values and initial conditions, the now pseudo-Goldstone boson 
can act as a dynamical dark energy component. In this new situation, the case /? = will again be peculiar, 
because it is the only case where dark energy is purely dynamical and has no constant contribution. 

Based on these reasons, in what follows, we will single out the case j3 — {) and study the associated 
phenomenology in more detail. 

B. Combining scale invariance and Unimodular Gravity 

We now want to add to the idea of SI the idea of Unimodular Gravity (UG) ]23-30] and apply it to the 
Higgs-Dilaton scenario. In UG one reduces the independent components of the metric g^^i, by one, imposing 
that the metric determinant g = dei{g^^) takes some flxed constant value. Conventionally one takes \g\ = 1, 
hence the name. Fixing the metric determinant to one is not a strong restriction, in the sense that the family 
of metrics satisfying this requirement can still describe all possible geometries. 

If we impose the unimodular constraint, the scale-invariant Lagrangian (5) becomes 

Csi^UG =\ (CxX' + 2a<^V) ^ + ^SM[A^o] - \r''^^,x^.X - V^'^. X) , (H) 
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where a hat on a quantity, like R, indicates that it depends on the unimodular metric g^j^, which satisfies 
detg^i, — —1. The potential l^(<p,x) is still given by (6). As a consequence of the unimodular constraint, 
the action associated to Csi-UG is no longer invariant under all diffeomorphisms (Diff), but only under 
transverse diff'eomorphisms (TDiff), i.e. coordinate transformations a;^ h- > x^ +^^(a;), with the condition 
d^^^ = 0. Just as in pure UG (not including non-gravitational fields) the equations of motion derived from 
the Lagrangian (11) contain an arbitrary integration constant Aq, which can be interpreted as an additional 
initial condition. It was shown in [1] that the classical solutions obtained from the Lagrangian (11) are 
equivalent to the solutions derived from the equivalent DifF invariant Lagrangian 

rSI-UG 1 1 

(exX' + 2a^V)^ + ^SM[A^o]-;^.9^''5^x9.x-V'(^,x)-Ao, (12) 



where Ao is the mentioned arbitrary constant. While in the original formulation (11) the dimensional 
constant Ag only appears in the equations of motion and thereby spontaneously breaks SI, in the equivalent 
Diff invariant formulation (12) the same constant appears as an explicit symmetry breaking in the action. 
Nevertheless, this constant should not be understood as a parameter in the action, but rather as an arbitrary 
initial condition. Given the equivalence of the two formulations, in order to study the phenomenology issued 
by (11), we will simply study the theory given by (12) for different values of Aq.^ 

We now turn our attention to the physical implications of the term proportional to Ag. In the first instance, 
let us only consider the gravitational and the scalar sector of the theory, i.e. 

^ \{^xX^+^hh')R-\{d^xf-\{d^hf-V{h,x)-^o. (13) 



y=(/ 2 ^'^'^ ^" ' 2' ^'^' 2 

where h is the Higgs field in the unitary gauge. In order to simplify the physical interpretation we define 
the Einstein-frame (E-frame) metric^ 

in terms of which the Lagrangian reads 

^ Ml^-\k-U{h,x), (15) 



V^g "2 2 

where if is a non-canonical but positive definite kinetic term (given below in (25)) and U{h^ x) is the E-frame 
potential given by 

Note that the E-frame potential gets singular at x = ^ = 0- The reason is that at this point the conformal 
transformation (14) is singular and the change to the E-frame is not allowed. Since for x — h = Q scale 
invariance is not broken, we will not be interested in the theory around this point. Let us discuss the shape 
of the E-frame potential and the classical ground states of the theory for a, A,^^,^/i > (cf. figure 1). If 
Aq — 0, in which case the theory (12) reduces to (5), the potential is minimal along the two valleys 

A A4x + oi.E,h 

They correspond to the infinitely degenerate family of classical ground states found in (7). As before, if 
/? = 0, the potential vanishes at its minimum, while a non-zero /? gives rise to a cosmological constant (9). 
In other words, spacetime in the classical ground state is Minkowskian, dS or AdS. These are the results we 



^ Of course, the choice to analyze the theory in the Diflf invariant rather than in the original formulation is purely a matter of 

convenience. 
* The Lagrangian in terms of the original variables is said to be written in the Jordan- frame ( J-frame) . 
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Figure 1: These plots show the shape of the E-frame potential U{h, x) (equation (18)) for Ao = 0, Aq > and Ao < 
respectively. 

have already discussed section II A. As soon as Aq 7^ the valleys get a tilt, which lifts the degeneracy of 
the classical ground states. For Aq < the valleys are tilted towards the origin. The true classical ground 
state for this case is the trivial one, x = ^ = 0. Hence, we discard this possibility. For Aq > the potential 
is tilted away from the origin, it is of the run-away type. In this case the theory has an asymptotic classical 
ground state, given by (17) with xo ~^ 00. Again, depending on the value of f3 this asymptotic solution 
corresponds to Minkowski, dS or AdS spacetime with curvature given by (9). 

We see that as a consequence of the non-minimal coupling between the scalar fields and gravity, the 
arbitrary integration constant Aq does not play the role of a cosmological constant (as it does in pure UG) 
but rather gives rise to a peculiar potential for the scalar fields. For Aq > the potential is of the run- 
away type. In the following sections we will see that such a potential can have an interesting cosmological 
interpretation. In fact, the evolution of the scalar fields along the valley can give rise to dynamical dark 
energy (quintessence). We will focus on the case /? = 0, where dark energy does not contain a constant 
contribution and is purely due to the term proportional to Aq (cf. arguments in section II A 3) . 

While the term proportional to Aq can play an important role in cosmology, its presence barely affects the 
particle physics phenomenology of the model. In fact, if the run-away potential is of the order of magnitude 
of the present dark energy density, the time evolution of the background scalar fields along the valley can be 
neglected on particle physics time scales. Also, the additional interactions between the Higgs field and the 
dilaton that are induced by this potential are negligibly small. 

C. Higgs-Dilaton Cosmology — The Qualitative Picture 

In this subsection we want to qualitatively describe the cosmological scenario issued by the model (12) (or 
equivalently (11)) presented in the previous two subsections (cf. [1]). We consider the theory in the Einstein 
frame (15) and focus on the case (3 = 0, for which the scalar field potential (16) reduces to (cf. figure 1) 



U{h,x) 



M* 



{^xx^+^hh^r 



'JH^ 



a . 



Ar 



(18) 



The scalar fields x ^^id h are now considered to be homogeneous background fields evolving in fiat 
Friedmann-Robertson- Walker (FRW) spacetime. Their evolution is affected by the non-canonical nature 
of the kinetic term. However, since the kinetic term is positive definite, in order to get a qualitative picture, 
it is enough to look at the features of the potential. In the absence of Aq, IJ has its minima along the two 
valleys \t? = f X^- The main effect of Aq ^ is to give a tilt to the valleys. As discussed in the previous 
subsection, Aq < is phenomenologically unviable. We will only consider the case Aq > 0, in which the 
valleys are tilted away from the origin. 

For an appropriate choice of parameters, the crude picture of the role of the cosmological scalar fields is the 
following: If the scalar fields start off far from the valleys, Aq can initially be neglected, and the scalar fields 
roll slowly towards one of the valleys. This roll-down can be responsible for cosmic inflation. As inflation is 
mainly driven by the Higgs field, this phase is much like in the case of the Higgs-Inflation model of [4]. 
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After the end of inflation preheating takes place. During this phase the scalar field dynamics is dominated 
by the field h and hence preheating in the present model is expected to be very similar as in the Higgs- 
Infiation model [31, 32] (see also [33]): The gauge bosons created at the minimum of the potential acquire a 
large mass while the Higgs field increases towards a maximal amplitude and starts to decay into all Standard 
Model leptons and quarks, rapidly depleting the occupation numbers of gauge bosons. The fraction of energy 
of the Higgs field that goes into SM particles is still very small compared to the energy contained in the 
oscillations, and therefore the non-perturbative decay is slow. As the Universe expands in a matter-like 
dominated stage with zero pressure, the amplitude of the Higgs field oscillations decreases. Eventually, this 
amplitude is small enough so that the gauge boson masses become too small to induce a quick decay of 
the gauge bosons. As a consequence, their occupation numbers start to grow very rapidly via parametric 
amplification. After about a hundred oscillations, the produced gauge bosons backreact on the Higgs field 
and the resonant production of particles stops. The Higgs field acquires a large mass via its interaction with 
the gauge condensate and preheating ends. From there on, the Higgs field as well as the gauge fields decay 
perturbatively until their energy is transferred to SM particles. 

The phase of preheating is followed by the usual radiation and matter dominated stages, during which the 
scalar fields are "frozen" at some point of the valley. Their energy density is now given by the Aq term and 
practically constant. As a consequence, it eventually comes to dominate over radiation and matter and hence 
provides a dark energy component. In other words, the scalar fields rolling slowly down the potential valley 
play the role of a thawing quintessence field [34-37[. In this late stage, the fields satisfy h{tY — f x(0^- 
On particle physics time-scales the time-variation of the background fields can be neglected. Perturbations 
around this almost constant symmetry breaking background can be interpreted as the SM particles plus an 
additional almost massless and almost decoupled particle, the dilaton. Note that as long as the background 
is constant, it is equivalent to quantize perturbations in the original (Jordan-) frame or in the Einstein- frame 
(cf. [38, 39[). 

In the following sections we present a detailed analysis of the infiationary phase (section III) and the 
dark-energy dominated phase (section IV). A detailed study of preheating in the present model is left for a 
further work. 

Allowing for /3 7^ in the potential would not affect the discussion of inflation. It will, however, have an 
effect on the dark energy phenomenology, on which we will comment in section IV. 

III. HIGGS-DILATON INFLATION 

As usual, it is assumed that during inflation all the energy of the system is contained in the inflaton flelds 
and in the gravitational fleld. Therefore, during this stage, the SM fields can be neglected. Let us rewrite 
the scalar-tensor part of (12) as 

-^ = ^-^R - \9^''5ai,d^rdA' - UW , (19) 

with a non-minimal coupling 

/(0) = ^ea</'"', (20) 

a 

and the potential 

[/(</)) - V{cb) +Ao = ^(^h'- ^x') ' + Ao , (21) 

including the SI breaking term Aq. As discussed in section HA, the parameter a is set to be very tiny 
a ~ 0(10""^"), in order to obtain the correct hierarchy between the electroweak and the Planck scale. Greek 
indices ^,1^, ... = 0, 1,2,3 denote spacetime coordinates while Latin indices are used to label the two real 
scalar fields present in the model: the dilaton field (p^ = x ^^d the Higgs field in the unitary gauge (j)'^ = h. 
The abstract notation in terms of 0* will in the following allow us to interpret the scalar fields as the 
coordinates of a two-dimensional sigma-model manifold. We will be able to write expressions and equations 
that are covariant under variable changes (f) i-> (p' {(j)). 
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Whenever the non-minimal coupling is non-zero^ /(</)) ^ 0, one can define the new metric 

~9^u = f^'gp. , (22) 

with O^ = Mp f{(f>) to reformulate the Lagrangian in the E-frame. Taking into account that the metric 
determinant and the Ricci scalar transform as^*^ 

^/^ = n-^y^ , (23) 

R = n^ (R + 6n\nn- Gg^^'df, In ft d^ In n\ , (24) 

one obtains 

-^ = ^R-Ik-U{^), (25) 

\/—Q 2 2 



where the kinetic term is given by 



K - labr''^^.r^,4>'' , (26) 



and "fab is a generally non-canonical and non-diagonal field space metric, which in terms of the variables 
{(j)^,(f>'^) = (x, /i) is given by 



02 \— ■ 2^ f^ 



7a.= ,^Ua.+ ;,M^^^ . (27) 



Unlike in the single-field case, the non-canonical kinetic term can not in general be recast in canonical form 
by redefining the scalar field variables. In fact, the field-space metric can be brought to canonical form 
by a local variable change if and only if its Riemann tensor identically vanishes. In the present case of a 
two-dimensional manifold, the Riemann tensor has only one independent component, and it is enough to 
compute the Ricci scalar Rj associated to the field space metric jab, 

For R^ to vanish globally, one would need to have ^^ — ^h- As we will see, this case is not allowed by 
phenomenology. The E-frame potential is defined as 

U{^)^V{<P) + Va,{<P), (29) 

where we have defined a scale-invariant and a scale-invariance breaking part as 

n<P)^^ and Va„(.^) = ^. (30) 

We can now write down the equations of motion derived from the E-frame Lagrangian (25). Einstein's 
equations are 

G^. = Jab (d^rd.<l>'' - \9^,r^dprdacl>'') +U9^.u, (31) 

where G^^ is the Einstein tensor computed from the metric g^^. The equations for the scalar fields are 

n<P'+r''Kbd^.rd,<P'^U-'\ (32) 



® For our choice of parameters, where ^x'^h > this is the case whenever the scalar fields are away from the origin (x, h) ^ (0, 0). 
^^ The action of the covariant d'Alembertian D on a scalar field s[x) is given by Qs = j-^ 8p (^yZ—gg'^'^dus). 
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where F^^ is the Christoffel symbol computed from the field-space metric ^ab, 

Tafc = 2^"'^ ('^'^'"'^ + ^d.^^'' ^ ^°-^^d) , (33) 

and where we use the notation JJ'''^ = ^'^'^U,d- Notice that equations (31) and (32) are covariant under 
redefinitions of the scalar field variables (ff i—)- 0'* (</)). 

We choose to do our analysis in the Einstein-, rather than in the Jordan frame. The reason for this choice 
is that in the literature predictions for measurable quantities are usually computed in the Einstein frame, 
where gravity has the standard GR form. At the classical level there is, apart from such practical arguments, 
nothing that would privilege one or the other frame. After all, the choice of the frame simply corresponds 
to a choice of variables. 



A. Exploiting scale invariance 

By construction, all terms in the Lagrangian (19), except the one proportional to Ao, are invariant under 
the scale transformations (4) . We will see that if the Ap-term is to be associated with dark energy, it must be 
negligibly small during inflation. The approximate scale-invariance of the theory will considerably simplify 
the analysis of the inflationary period. 

1. The Noether current of scale invariance 

Let us start by computing the Noether current associated to the scale transformations (4), which for an 
inflnitesimal value of the parameter a become 

9t,u ^ .g^i. + crAg^,, , 0* 1-^ (/.* + crA0* . (34) 

The explicit expressions for Ag^^ and Ai/i* depend on the choice of the fleld variables. For the original 
variables one has Ag^^ — —2g^^, Ax = x ^^^d A/i = h. The associated current is given by (see e.g. [40]) 

dC dC 

V-gJ'^ = J^ -,^9c.p + -WT^r^A^' ■ (35) 



and satisfles 



1 9[Ao7= ^] 
^ dg,iu 



D^J^ = - ^ ^ "^ ^^ Aff^. = 4Ao , (36) 



where -D^ denotes the covariant derivative constructed with the metric g^i^. 

In the E-frame, scale transformations do not act on the metric, Ag^i^ = 0, and are simply given by 

0* h^ (/)' + ctA(/)^ . (37) 

In this case, the expression for the current is 

while the conservation law becomes 

where the covariant derivative _D^ is constructed with the metric g^^. Whenever Aq vanishes, scale invariance 
becomes exact and the associated current J^ , respectively J^ is conserved. 
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2. New variables 

The approximate conservation law can lead us to a very convenient choice for the scalar field variables 
in the E-franie formulation. In fact, one can always choose a set of variables (</)' ,<p' ) — (p, 6) such that a 
scale transformation only acts on a single variable p. Moreover, p can always be defined such that the scale 
transformation acts on it like a shift p i-^ p + aMp. This transformation is a symmetry of the Lagrangian 
if Ap = 0. This means that p can appear in the Lagrangian only through Vaj, and through derivatives in 
the kinetic term. Moreover, the variable 9 can be defined such that the field-space metric 7^^, which is 
independent of p, is diagonal. If the Lagrangian (25) is expressed in terms of variables that satisfy these 
requirements, the current J^ takes the form 

> - Mpr''ipp{e)d,P ■ (40) 

Now, in order to find the relation between the variables (0' , 0' ) = (p, 9) and the original variables {(}?-, (fp') = 
(x, h), let us express the current J^ in terms of the original variables. 

Comparing this expression to (40), it is clear that p must be a function of the combination of the fields x ^^nd 
h acted upon by the partial derivative, i.e. p = p[{l + 6^^) X'^ + (1 + Q^h) h"^] ■ For the scale transformation 
to correspond to a shift of p by aMp, this function has to be chosen as 

Mp^ f{l + 6^^)x^ + 0^ + 6^h)h^ 
Pp 

The variable 9, as it does not transform under scale transformations, has to be a function of the ratio 



"-T^'- r "" Mi 1^ c^' 



between h and %, i.e. 9 — 9 



There is some freedom in the choice of this function. One can notice that 



h 
_X_ 

the argument of the logarithm in (42) corresponds to the radius of an ellipse in the (x, /i)-plane. We will 
therefore define 9 as the angular coordinate of the ellipse, i.e. 



Let us note that in terms of the variables {p, 9), since one has Ap ~ Mp and A9 — 0, the current conservation 
law (39) corresponds to the equation of motion for p. Further, one can see from (39) that the dependence 
of Vaj, on p is such that Vaq oc exp (— 4p/Mp). 

In terms of the new variables, the E- frame kinetic term K and the potential U = V + V^g (cf. (25)) are 
given by 



and 



£.h J sin^ 9 + ^ cos2 9 ^x cos2 9 (tan^ 9 + <;y 



XMt /sin20-fj±§^cos2fl\' 



^(^)-^ sin^;;::tsv i ' (^^) 



'i + ee 



2 „-4p/Mp 

T/A„(p,e)=Ao(^^^) — — —,, (46) 

(sm"^ 9 + <; cos2 9) 



where we have defined the parameters 



Cx _(i + 660ex 



^^f' ^^ : :r / ■ (47) 

^h (1 + 6^x) ^h 
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We will see in section III D that for a successful description of inflation the parameters have to be such that 
^^ ^ 0(10""^) and ^/j ~ 0(10"'') and hence /i ^ 1. In this case, one can neglect ^ in the kinetic term (44). 
In this approximation, the action can be further simplified by introducing the variables 

p^l^^p and e= ^- tanh"^ Ul - <; cos 9) , (48) 

a ^ ' 

with the parameters 




^"^^ '^ 7=Wr^^- (49) 



^ V 1 + eCx 

In terms of these variables and for /x <C 1, the kinetic term takes the simple form 

k ^ e^^'^'^\d~pf + {d~Bf , (50) 

with 

h(d) = - In r^cosh^ [oBlM^ , 
which has been studied in the literature previously [41-43]. The potentials are given by 



V{B) - ,^;;:-^^ {\ - ^cosh^ (a^VMp) - '^\±^^x^\?{aBlMp)\ , (51) 






Va„ (p, ^~) - ^ <r' cos\,\ae/Mp)e-^'^p/'''- . (52) 

5. Departure from scale invariance 

Let us now look at the E- frame equations of motion (31) and (32) in order to see in which region of field 
space the effect of a non-zero Aq will be important. Ao enters the equations through U — V + V\g and 
through U''^ — V''^ + Vj^ . We therefore define the two new parameters 

vi^^, (53) 



V2 = I — , (54) 



that characterize the departure from scale invariance. The parameter vi compares the importance of the scale 
invariance breaking part of the potential Vao to the scale-invariant part V in the Einstein equations (31). In 
the region of field space where vi ^ 1, Vaq can be neglected in Einstein's equations. The parameter V2 pro- 
vides a coordinate invariant measure of the importance of Vao'" compared to V'"" in the scalar field equations 
(32). In fact, locally one can always choose a coordinate system in which |Vao'"'^/V^'"'^| == I^Ao'^/^'^l — ^2- 
Hence, in the region of field space where V2 <C 1, Vaq'"' can be neglected in the scalar field equations. In 
the region where both vi ^ 1 and ^2 ^ 1 hold, the effect of Aq is negligible and the equations of motion 
become practically scale-invariant. We will refer to this region as the scale-invariant region. It will turn out 
that for phenomenologically viable values of the parameters, the whole period of observable infiation takes 
place in this region. 

B. Evolution of the homogeneous background 

Let us now consider homogeneous scalar fields 0* = 0*(O i'^ spatially fiat FLRW space-time characterized 
by the line element 

ds^ = gf.^dx^'dx" = -dt^ + a^{t)dx^ . (55) 
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Before writing down the equations of motion, we introduce some notation for vectors lying in the tangent 
and cotangent bundles of the field manifold, that will allow us to write many of the upcoming expressions 
in a very compact way. The notation corresponds to the one of [44]. We denote vectors in boldface, i.e. 
A = (A^jA^). The inner product of two vectors A and B is given by 

AB = A^B = "i^jA'B^ , (56) 

and the norm of a vector A is 

\A\ = VA-A, (57) 

where a dagger ^^ on a naturally contravariant or covariant vector denotes its dual, e.g. 0^ = {jij(p-') and 
V^ = (7*-'Vj). Here and in the following we use Vj to denote the covariant derivative constructed from the 
field-space metric (7^^). 

For homogeneous fields in flat FLRW space-time, equations (31) and (32) reduce to the Friedmann equa- 
tions and the equations of motion for the scalar flelds, 



^-3M| Ul'^l +^1' ^''^ 



2i? + 3i7^^-^(l|0r-C/), (59) 

^+3H0 = -V^t/, (60) 

at 



where a dot stands for a derivative with respect to t, H = a/ a and the action of Z) on a contravariant vector 



X' is deflned as DX^ — dX'^ + F^f^X^d^''. To this we can add the equation for the current (39), which for 



homogeneous fields reduces to 

|(a30.A0)=4VA„, (61) 

where A0 = (A0^, A(/)^). This relation is of course not independent of the equations of motion. However, it 
will prove useful in the following. 

For the discussion of inflation, it is helpful to change the time parameter from t to the e-fold time parameter 
N — Ina(i). The field equations can then be written as 

(62) 
(63) 

cj)' = -MlV'^ In U , (64) 

3-i|0f/M2 

where a prime denotes a derivative with respect to the e-fold parameter N . 

1. Slow-roll inflation and background trajectories 

In the present model, inflation can occur due to a phase of slow-roll of the scalar flelds over the flat part 
of the potential towards one of the potential valleys (cf. flgure 1). Let us deflne the slow-roll parameter e 



H^ 


U 


■AMI -\(t>'^ 


H' 


10'^ 


~H 


2 Ml ' 




D^' 




dN. + cb' 



and the slow-roll vector rj as [44, 45]^ 



H' _ 1 |0f 
~H ^ 2 Af J, ' 
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(65) 



-^/l^'l' (66) 

in terms of which equations (62) and (64) read 

-^|^'|+0' = -M|,V^ln[7. (68) 

O 6 

The exact condition for inflation, i.e. for a > 0, is given by e < 1. 

The slow-roll regime is characterized by the fact that the equations (67) and (68) are well-approximated 
by the "slow-roll equations" 

U 



H^ = -^ , (69) 

3M|, ' ^ ' 

(j)' = -MlvHnU . (70) 
The conditions for the validity of the slow-roll approximation are^^ 

e < 1 and 77 = [r?] < 1 . (71) 

Still following ref. [44] we introduce the kinematical orthonormal basis vectors en = t^, pointing in the 

|0'| 

direction of the field trajectory, and ej^, pointing in the direction of (/ — e||en)?7, where I is the 2x2 unit 
matrix. This allows us to write the second slow-roll parameter as 

V = Jril + vl , (72) 



where the speed-up rate m and the turn-rate ryj^ are defined as 

77|| =e|| -7) = -^^ {\(j)'\+Mle\\ •Vint/') , (73) 

r?_L = e_L • rj = -Ml -p^ e_L • V In t/ . (74) 

|0'| 

In the case of one-field inflation, ryj^ is equal to zero and 77 = \'ri\\\. In the upcoming section we will see that for 
the Higgs-Dilaton model, as long as the flelds are in the scale-invariant region, ryj^ goes to zero very quickly, 
and rj can be computed like in a single fleld model. 

Within the slow-roll approximation, i.e. making use of equations (69) and (70), one can compute approx- 
imations in terms of the potential for the slow-roll parameters e and 77 

e(SR) = lM2|Vlnt/|', (75) 

^(SR) ^ /(e(SR))tM2ef«') , (76) 



^^ The definition of the vector rj used here differs from the definition in [44] by the factor \4>'\- 

^^ Note that e = ei and |»?| = 2 kzl, where ei = -fr and €2 = —j^ are the multifield generalizations of the standard first two 
horizon-flow parameters defined in [46]. 
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such as for the speed-up rate r/n and the turn-rate rj± 

,?i^«')=-(epVMep, (77) 

^(SR) ^ _(g(SR)^tMeV''^\ (78) 

where the matrix M is defined as 

M = M|.V'I' Vint/ (79) 

and the kineniatical unit vectors in the slow-roll approximation are given by 

,SR)^_Vnnt/ ^^^ _^M+£^)^ 

Hence, instead of the exact slow-roll conditions (71) one can use the approximate slow-roll conditions 

e(SR)«i and r^^^R) « l , (81) 

which should be understood as consistency conditions for the slow-roll approximation. Once the system is in 
the slow-roll regime, i.e. the exact conditions (71) are satisfied, the approximate conditions (81) guarantee 
that the system remains in the slow-roll regime and describes a phase of infiation. We will approximate the 
time where infiation ends, as the moment where e^^^' — 1. 

Let us now discuss the regions in the (x, ft,)-plane for which the approximate slow-roll conditions hold (cf. 
Fig. 2) . The slow-roll region extends to infinity along the potential valleys if ^^ < ^ . As will be shown in 
section IV, only if this condition holds, the scalar fields can constitute a dark energy component in the late 
stage. Further, during infiation it is safe to neglect the term in the potential proportional to a. In fact, 
for a = the potential possesses only one valley which goes along the x-axis. For a 7^ this valley splits 
into two valleys that lie at the angles 6 = zt arctan(a) with respect to the x-axis. For a ^tC 1 these angles 
are very small. We will see that infiation in our model occurs far from these valleys where the effect of a 
non-zero a is irrelevant. Hence, we will put a = for the rest of this section. The plot of the slow-roll region 
for £,x < -^ and a = is presented in Fig. 2. 

Next, we want to analyze the different trajectories the fields can take if the initial conditions are chosen in 
the slow-roll region. We will only consider trajectories starting in the first quadrant x? ^ > 0. Trajectories 
starting in other quadrants are exactly analog. The shape of the potential (29) makes that all trajectories 
tend to approach one of the potential valleys. There are trajectories (type a) that on their way to the valley 
never leave the slow-roll region. Numerical computations show that such trajectories undergo only very few 
slow oscillations before asymptotically approaching the valley. One can not expect a successful reheating 
period from this type of behaviour [31, 32]. The good trajectories (type b) are those that at some point 
leave the slow-roll region. After the exit of the slow-roll region, which at the same time marks the end of 
infiation, these trajectories undergo a fast-roll towards the valley and therefore oscillate strongly around its 
minimum. Typical examples for both types of trajectories are given in Figure 2. 

Looking at Fig. 2 we observe that all good trajectories (type b, blue line) go through the scale-invariant 
region before leaving the slow-roll region. Therefore, for these trajectories, the end of infiation always takes 
place within the scale-invariant region. We will see in section IV B that requiring the scalar fields to act as 
a dark-energy component in the late phase will give a bound on the initial conditions. Qualitatively, this 
bound tells that the scalar fields during infiation have to be very far from the origin. Therefore, not only the 
end but the whole period of observable infiation (i.e. the final ^ 60 e-folds) takes place in the scale-invariant 
region. This fact considerably simplifies the analysis. In particular, during infiation the scale-invariance 
breaking part Vaq of the potential (29) can be neglected, i.e. [/ ~ V. As a consequence, if one uses the set of 
variables (p, 9) introduced in section HI A 2, the potential only depends on 6. Consequently, also the slow-roll 
parameters e^^^' and r/'^^' are functions of 6 only. Further, one finds that the lowest order approximation 
of the turn rate r]j_ vanishes, i.e. rj]^ = 0. 

In the scale-invariant region and written in terms of the variables {p,0), the slow-roll equations for the 
scalar fields (70) read 

p'-O, (82) 

e' = — ^%^ cot e (1 + ;^^ — ^] . (83) 
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Figure 2: The blue region is the slow-roll region for ^^ <§:. 1, ^h 2> 1 and a — given by e'^^' < 1. The inclusion of 
the second slow-roll condition t]^ ' < 1 does not change the essential properties of this region. The general features 
of the slow-roll region are the same whenever Cx < | ^^'^ C^ > i • For ^;,f < | and ^^ < | the central fast-rolling 
region vanishes. For ^^ > | the slow-roll region does not extend to infinity along the x~^^^^ in which case the scalar 
fields can not act as dark energy in the late stage of evolution. The shaded region corresponds to the scale-invariant 
region delimited by ui < 1 and V2 < 1. This is the region where the influence of Ao 7^ is small. The presence of 
the slow roll-region along the x-axis such as the central fast-roll region are effects of Ao > 0. For Aq = the slow-roll 
region is simply given by the triangles delimited by the two diagonal lines. Note that in this case the units of the axis 
have to be chosen differently. The red line represents a trajectory of type a, never leaving the SR-region. The blue 
line is a trajectory of the type b, which leaves the SR-region and oscillates strongly before rolling down the valley. 
These trajectories were found by numerically solving the exact equations (60). 



From the first equation one finds that tlie background trajectories in the scale-invariant slow-roll region 
correspond to a constant value for p, p ~ po, respectively to ellipses in the (x, /i)-plane described by 



(1 + e^x) x' + (1 + 6^) /i' = AfJ, e^^o/A^^ 



(84) 



The second scalar field equation (83), as a consequence of scale invariance, does not depend on po- This 
equation can be integrated in order to get the number of e-folds before the end of inflation as a function of 
the angle 9, 



N{9, 0end) = 



1 



In 



cos Wend 



+ 3^;, In 



S,^ COS^ ^cnd + ^h sin 6'ond + Q£.xOi 



cos« J - \ ^^ cos2 6* -f- C/i sin^ 6' + 6CxCh 
where ^end is the value of at the end of inflation, which can be found from the condition 

8^^(1 + 6^) COt^^end 



e^^^n^cnd) 



1 + 6Cx ^y C0S2 eiend + ^h sin 6'end 



1 



(85) 



(86) 



After inserting values for 1^^, S^h and requiring a minimal number of inflationary e-folds A^(0min, ^cnd) — -^min, 
equations (85) and (86) can be solved to obtain a lower bound ^min < ^initial on the initial conditions for 
inflation. In what follows, we will derive bounds on the parameters ^^ and ^^, which are related to the 
spectra of primordial perturbations. 
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2. One-field attractor in the scale-invariant region 

In the previous section we have seen that if the slow-roll conditions hold and if the system is in the scale- 
invariant region, i.e. vi ^ 1 and V2 ^ 1, the trajectories are given by p = po- In other words, there exists a 
set of variables in term of which only one of the two fields evolves during infiation. 

In this section we are going to show that in the scale-invariant region p' = is an attractor independently 
of slow-roll and also after inflation. We will further show that during inflation in the scale-invariant region, 
the turn-rate rjx^ goes to zero very rapidly, also beyond the lowest order slow-roll approximation. This fact 
will have important consequences for the evolution of perturbations of the scalar fields. 

Let us look at the equation for the scale-current (61). In the scale-invariant region, the term on the 
righ-hand side can be neglected. In terms of the variables (p, 9) and using the e-fold time parameter, one 
obtains 

P' - p^e-'^ , (87) 

where the constant depends on initial conditions, jpp can be read from (44) and represents a bounded 
function. During infiation the factor H~^ is generally nearly constant but grows at most like e^ at the end 
of infiation. For a matter domination like stage one obtains H^^ oc es" and for a radiation domination like 
stage H^^ ex 63^. One can read from equations (62) and (63) that whenever the potential is positive, which 
is always the case in the scale-invariant region of our model, H^^ grows more slowly than e^^. We conclude 
that as soon as the system can be approximated by the scale-invariant equations, p' = is an attractor. 
Let us now turn our attention to the turn-rate ri±. In the scale-invariant region, equation (74) reduces to 



r;_L = -Ml ^^ e_L • VlnF = -MJ,\VlnV\ ^. — ^ e_L 



e 



(vy) 



where e/^y% is a unit vector pointing in the direction of W. From the current conservation law (61), still 
neglecting Vaq and noticing that A0 is perpendicular to VV^, we find 



est. _3JY 



where e^^y-, is a unit vector perpendicular to e^^y-, . At every point in field-space one can find an orthogonal 
matrix O such that ej_ = Oen and e,^y-, ~ Oe,.^y\ . Hence, one has 



ej^ • e 



(vv-) 



= (Oe_L)- Oe.^^J = fiwe,^y. . (90) 



(vy)j 



Inserting this relation and (89) in (88) and making use of (65) and (75), we find the following result for the 
turn rate 

where we have redefined the constant related to initial conditions. During slow-roll infiation, e is growing, 
but remains smaller than unity, H is decreasing very slowly and e(^^) is small. |A0| is a bounded function, 
which in terms of the variables (p, 0) is given by |A</>| — Mp^Ffpp. We conclude that during slow-roll infiation 
and for approximate scale invariance, the turn rate r\\^ goes to zero exponentially fast. Further, one can 
show that if scale invariance is slightly violated, 77j^ is proportional to v^. 

C. Linear perturbations 

The theory of cosmological perturbations as stemming from quantum fiuctuations during infiation was 
developed in [47-51] (see also [52] and references therein). Including scalar and tensor perturbations and 
fixing the Newtonian transverse traceless gauge, the line element can be written as 

dip- = - (1 + 2$) dl^ + aitf ((1 - 2*) 8ij + h^^) dx'dx^ . (92) 
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$ and ^ are the Bardeen potentials [53]. For comparison with CMB observations, we will be interested in 
the power spectrum of the comoving curvature perturbation, which is defined as [53, 54[ 

C ^ * - ^ f ^ + H*j . (93) 

Through the perturbed Einstein equations, ( is related to the linear perturbations of the scalar fields, S(f>, 
like (see e.g. [55[) 

I0'l 

where the quantity in parenthesis is the multifield version of the Mukhanov-Sasaki variable [56, 57[ The 
evolution equation for ( is given by [42-44[ 

^'-M7)W^*-'^^M' ^''^ 

where A ~ S^^didj. The quantity 6(j)j_ = S(f> • ej_ is the component of the field perturbations perpendicular 
to the background field trajectory sometimes called relative isocurvature (or entropy) perturbation (see e.g. 
[58[). In the long wavelength limit, k <^ aH, the first term in (95) can be neglected and the evolution 
equation becomes 

This is a well-known result, showing that for multi-field inflation, C is not in general conserved outside 
the Hubble horizon [42, 43[. There are two cases in which the source term on the right-hand side of the 
equation vanishes. One is if the perturbation vector d(j) is tangent to the field trajectory, i.e. 54>± — 0. 
This corresponds to the complete absence of relative isocurvature perturbations during inflation and is 
not satisfied in our scenario. The second possibility is to have the potential gradient VC/ parallel to the 
background field trajectory, resulting in the vanishing of the turn rate ri± = 0. In our model, as we have 
seen in the previous sections, inflation takes place in the scale-invariant region, where ?7_l = is an attractor. 
Hence, as a consequence of scale-invariance, the comoving curvature perturbation is practically conserved 
outside the horizon, just like in single-field infiation.^'^ As mentioned before, corrections due to deviations 
from scale invariance are suppressed by the parameter f 2 . This parameter being extremely small (cf. section 
IV B), we will neglect such corrections. 

We will from now on use the fact that infiation takes place in the scale-invariant region and suppose that 
initial conditions are such that the attractor r]± = has been reached before the observable scales cross the 
horizon. In that case, ( is constant outside the Hubble horizon and the result for the primordial spectrum 
of ( is like in the case of one-field infiation. To lowest order in the slow-roll parameters it can be expressed 
as [52, 59, 60[ 

1 /H*V 1 V* 



^^(^^ - 2M2e* [2tt) - 247r2e(SR)* Mf, ' ^^^^ 

where quantities with an asterisk are evaluated at the moment of horizon crossing, i.e. when aH — k. The 
scalar spectral index is given by 

ns{k) - 1 - '^ - -2(.* + ryP ^ -2(e(^«)* + ,f ^)*) , (98) 
while the running of the spectral index can be expressed as [61[ 

ac(fc) . ^ . -Ae*,l - 2,IC c. -4.(-)*,i-)* - 2 (,^r) ^'"^ , (99) 



^^ Using the results of [44], one finds that 5(^x/l</''l is of the same order of magnitude as <^ and does not grow considerably 
during infiation. 
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where ^ is the third Hubble-flow parameter ^ = -jj^- In the slow-roll approximation the combination r/||^ 
can be expressed in terms of the potential as [44] 

with the matrix X = Mp'V In t/V^ Af , where U can be replaced by V due to scale invariance.^^ 

Now, even if the relative entropy perturbations 54)±^ do not affect the evolution of C,, since they are 
completely decoupled, they can in general be present at the end of inflation. The total entropy perturbations, 
however, are given by 5 ex ri^54)^/\4''\ [58]. And, as we have shown above, scale invariance during inflation 
leads to ry_L ~ and hence to a strong suppression of S. As has been shown in Ref. ]62] (see also [63[ 
for a general perspective), the large scale suppression of entropy perturbations during inflation avoids the 
resonant growth of these fluctuations also during (p)reheating. Working with this assumption, we will be 
able to relate the primordial spectra to CMB observations. 

The primordial spectrum of the tensor perturbations is given, to the lowest order in the slow-roll approx- 
imation, by [52, 59[^^ 

^.W-xftUt -:^77r> (100) 



Ml \2Tr J 37r2 M^ 



which results in a tensorial spectral index 



The ratio of the tensor and the scalar spectra to first order in slow-roll is then given by 

r= 5^~16e*~16e(s^)*, (102) 

and we have the consistency condition like in one-field infiation 

r = -8ng , (103) 

valid to the lowest non-trivial order in the slow-roll approximation. 

D. CMB constraints on parameters and predictions of the model 

In this section we are going to explicitly compute the primordial spectra and confront them with CMB 
observations. As discussed in the previous sections, the whole period of observable infiation takes place in 
the scale-invariant region. We can therefore use equations (85) and (86) for the background. Moreover, due 
to scale-invariance, C is conserved for large wavelengths during infiation. Further, we make the assumption 
that after infiation entropy perturbations die away before having an observable effect. This allows us to 
directly compare the primordial spectra (97) and (100) to observations of the CMB. 



Let us mention again that, since inflation takes place in the scale-invariant region, one can choose variables (p, 6) for which 
the approximate slow-roll parameters e'^^', 'q\ and the combination ('?||C) depend only on d. In models of one-field 

inflation with a canonical kinetic term an alternative common definition of approximate slow-roll parameters is given by 
es = 2^P ( TT ) ' Vs = ^p^r ^^'^ 5« — ^p jj-2 ■ In the case of one or several fields with non-canonical kinetic terms, 
these definitions generalize to es = j-^^pl^lnf^l > Vs = {^« V^sSu and gs = (ej, )^Xse\ where the matrices 
Ms and Xs are defined as Ms = M|, jjVtVC/ and Xs = M|,yj2 Vt/VtV+Vf/ and the unit vectors are given in (80). These 
parameters are related to the approximate slow-roll parameters used in the present work as t^°^' = es, 77^ ' = 2es — r]s and 

(,,||5)(SR)=8e2-6es,,s-h6- 

This result is based on the slow-roll approximation and involves no further assumptions. 
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We will show in the following, that the running of the scalar spectral index and the amplitude of tensor 
perturbations are related and very small, cf. (119). Therefore, we can consider those observational bounds 
(WMAP7 + BAO + Hq) for the scalar tilt and the amplitude of the scalar power spectrum, which are based 
on the standard ACDM model and the assumptions that the primordial spectrum obeys a power-law and 
that tensor modes can be neglected (see Ref. [64]) 

Vcika) = (2.43 ± 0.27) X 10"^ (104) 

n,(fco) = 0.968 ±0.036, (105) 

where ko/ao = 0.002 Mpc^^ and the indicated errors correspond to the 99% confidence levels. 

Let us start by computing the spectral quantities Vi^iko), ns{ko), a(fco) ^^nd r(fco) evaluated at the pivot 
scale fco, in terms of the parameters ^^, ^/j and A. This is done in four steps: 

i. Equation (86) is solved for 9end = d end{£,x^ ih) ■ 

ii. Oend is inserted into (85) from which one determines 6* = 9*(^^,^h, N*). 

iii. Expressions (97), (98), (99) and (103) are evaluated at 9* to find the spectral quantities as functions 
ofex,a, AandTV*. 

iv. N* , the number of e- folds between the moment where fcg exits the horizon and the end of infiation, is 
expressed as a function of the parameters ^^, ^h and A. 

In order to determine N* (step (iv)) we need to know the post- inflationary evolution of the universe including 
the details of the reheating process. If there are uncertainties related to the post-inflationary history, these 
can be accounted for by varying the value of N*. One can compute N* approximately, by making a few 
assumptions about the post-inflationary evolution. First, during the reheating phase the scale factor is 
expected to evolve like in a matter dominated universe. The reason is that during this stage the present 
model behaves much like the Higgs-Inflation model (cf. [31, 32[). In this sense matter-like scaling of the 
universe during reheating is not really an assumption but rather a property of the considered model. Next, 
we make the usual assumptions that reheating is followed by the standard radiation and matter dominated 
stages. Further assuming that the transitions between the different phases are instantaneous, one can derive 
the following relation (cf. [59[) 

Here, ao, Hq, q'^ and I^q stand for the current values of the scale factor, the Hubble parameter, the critical 
density and the abundance of radiation, respectively, g^h denotes the radiation energy density at the end of 
reheating, i.e. at the onset of the hot big bang. After inserting the observational value figft,^ ~ 4.2 • 10^^ 
(for T^ ~ 2.7ZK [65[), where h is the dimensionless Hubble parameter, the above formula can be written as 

iV*.59-ln^^^^-ln^r^+ln(.^^^) _ . 1, , ^^^ , . (107) 

0.002ao y(6i*)V4 > --'- ^ ' o I . I v ; 

Notice that the dependence on Hq has cancelled out^^. A detailed determination of Qrh goes beyond the 
scope of this work and is postponed for a future publication. We can, however, consider two limiting cases. 
An upper limit on g^.^ is simply given by 

e™- = F(0e„d) , (108) 

corresponding to instantaneous reheating at the end of inflation. A lower limit can be found in the same 
way as for the Higgs-Inflation model in [31, 32[. Namely, one can look for the value of 0, below which the 
particle interactions become those of the Standard Model, apart from the suppressed interactions discussed 





In the analog formula of [59] this fact remains somewhat hidden. 
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Figure 3: This plot shows the parameter regions for which the amplitude P^{ko) and the tilt ns{ko) of the scalar 
spectrum lie in the observationally allowed region (WMAP7 + BAO + Hg at 99% confidence level), for A = 1. (The 
variation of the result induced by variation of A in the interval 0.1 < A < 1 is negligible.) The red region is obtained 
for Qrh ~ (ffh^ (instantaneous reheating), while the blue region corresponds to Qrh = Q™"' (long reheating). The 
fact that the bands are cut on the right comes from the constraint on the scalar tilt ns(fco), cf. (116), while the 
band-shape is due to the constraint on the amplitude Pf (fco), cf. (115). 

in section II A 1 , and therefore guarantee immediate reheating. Inspecting the kinetic term (44) and the 
potential (45), we find that this happens as soon as t&v? 9 < min{ij,, 1).^^ We will therefore set the lower 
limit 



mm 
Qrh 



Trfnmm 



with 



tan^ 9, 



7mn 
h 



min{fi, 1) . 



(109) 



The value of A to be used when computing the spectral parameters corresponds to the Higgs self-coupling 
evaluated at the scale of inflation [31, 32]. It contains the uncertainty related to the Higgs mass mj^. We 
expect the running of A to be similar as in the case of the Higgs-Inflation model [66-73]. As was shown 
in [31, 32[, for nijj ~ 130 — 180 GeV, the coupling A evaluated at the scale of inflation lies in the range 
A -0.1-1. 

Step (ii), i.e. the exact inversion of (85) can not be done analytically. Hence, before computing approximate 
analytical results, we execute the four steps numerically. This allows us to plot the region in parameter space, 
for which V({kQ) and ns{ko) lie within the observational constraints (104) and (105), cf. figure 3. It turns 
out that in the observationally allowed range and for A ~ 0.1 — 1, the spectral quantites depend on £^h and 
A almost only through the combination ^^/vA. This fact will become explicit in the approximate analytical 
results to be derived below. The red region in figure 3 is obtained under the assumption of instantaneous 
reheating (grh = P^/f ^) while the blue region corresponds to the case of long reheating, i.e. grh 
obtain the bounds 



e;7,"'. We 



< Cx ^ 0.008 , 
43000 < 4^ < 85000 , 



for n 1, — „irim 
lOl grh — Qrh i 



(110) 



Note that the value of f hes always between 1 and the value of /i. 
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Figure 4: The spectral tilt as a function of the non- minimal coupling parameter ^^ . The other parameters are set to 
^h = 65000 and A = 1. Note, however, that changing the ratio 5h/\/A in the observationally allowed range affects the 
result only by a negligible amount. The solid curves correspond to the numerical results, the blue one is obtained 
for Qrh = Q^h" (long reheating) and the red one for Qrh = Q^it^ (instantaneous reheating). The blue and red dashed 
curves are obtained from the analytical approximation (116) for iV^j„ and N^^^. The black dashed curve represents 
the asymptotic solution (133), which is a good approximation if j^ < ^^ ^ 1 • The horizontal line and the shaded 
regions correspond to the observational mean value and the Icr and 3a confidence intervals, cf. (105). 



and 



< Cx < 0-008, 
44500 < ^ < 92000 



for Qrh = Q: 



max 

rh 



(HI) 



In this region of parameter space, the quantities ^^(fco), Q;(fco) and ^(fco) vary only with ^^. The numerical 
results for these quantities are given in figures 4 and 5, respectively. 

We have found that the parameters need to satisfy ^^ <C 1 and ^/i ^ 1. With this knowledge, we again 
carry out the above four steps and derive approximate analytical results. From (86), we obtain 



^end 



= 2 X 31- 



l + O 



X' 



(112) 



In order to approximately solve (85) for 
inversion then gives 



we can neglect the second term on the right-hand side. The 



arccos cos 



{Oend)e 



^tixN- 



(113) 



Here the sign for approximate equality "~" refers to the approximation made when inverting (85). This 
approximation will constitute the main source of error in the approximate expressions for the special quan- 
tities. One can get a more accurate approximation by reinserting the first approximation into the right-hand 
side of (85) in order to compute the second order approximation of an iterative solution. However, as the 
expressions get considerably more complicated, we stick to the first order approximation, which already 
comes very close to the numerical results. Inserting O^nd from (112) into (113), one obtains 



arccos e 



-ii.N'\ ^^0 



Cx. 



1 



(114) 
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Figure 5: The running of the scalar spectral tilt and the tensor-to-scalar ratio (right) as a function of the coupling 
£,-^. The other parameters are set to ^h = 65000 and A = 1. Note, however, that changing the ratio (,h/V^ in the 
observationally allowed range affects the result only by a negligible amount. Blue solid curves show the numerical 
results for Qth ~ qTh" (long reheating), while the red solid curves are the numerical results for the case Qrh = f?™^ 
(instantaneous reheating). The dashed curves are obtained from the approximate expressions (118) and (117), 
inserting iV;;^^ and N^^^ 



We can now evaluate the spectral parameters at the approximate value for 6*. Inserting (114) into (97), 
(98), (99) and (102) and recalling (103) we obtain^** 



Pdko) 



Asinh^(4^^iV*) 



l + O 



ex. 



n, (fco) - 1 ~ -8^x coth i4:^^N*) l + O 



■^X' e ' 



ac(fco) ~ -32^^ sinh"" (4Cx^*) U + C 



r(fco) = -8ng(fco) ~ 192^^ sinh"^ (4Cx^*) U + C 
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One can see that in this approximation, a(;{kQ), r(fco) and Ug^ko) are related as 



1 4 

-g'"(fco) = ^ng{ko) 



(115) 
(116) 
(117) 
(118) 

(119) 



which can be understood as an approximate consistency condition for the Higgs-Dilaton model. 

In order to find N* in terms of the parameters of the theory, we insert (114) into (107), from which we 
derive the approximate results^^ 



N* 



N* 

max 



64.3- ^lnA-^lni|)( l + O 



^X' ,- 1 



1 1 



64.5 - ^ In ^ ) ( 1 + O 



?X' e ' 



1 1 



(120) 
(121) 



^® Note that these approximate results can equivalently be derived from the approximate action given by (50) and (51) with 

a = 0. 
^® The numerical factor is given to the first decimal. 
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Figure 6: The number of efolds A''* between horizon crossing of the scale fco and the end of inflation as a function of 
ih/V^, for £,^ — 0.001. Note that changing ^-^ in the observationally allowed range affects the result by a negligible 
amount. The solid curves correspond to the numerical results, the blue ones correspond to N^i„ (long reheating) 
and the red ones to N^^^ (instantaneous reheating). The upper one of the solid blue curves is obtained for A = 1 and 
the lower one for A = 0.1, showing the slight dependence of N*„i„ on A. The dependence of N^^x on A is negligible. 
The blue and red dashed curves are obtained from the analytical approximations(120) and (121). 
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where the subscripts "min" and "max" stand for the cases p^ft, = fi"™" ^^nd Qrh — Q^h^ respectively, 
approximate results together with the numerical results for N*^^^^ and N^^^ are plotted in figure 6. 

Equations (115)-(118) together with (120) and (121) constitute our approximate analytical results for the 
spectral parameters in terms of the parameters ^^, ^/i and A, in the two limiting cases of instantaneous and 
late reheating. One can see that the results depend on ^h and A mainly through the combination ^/i/vA. 
Independent variation of A affects the spectral parameters only through N^^^^. Further, we remark that for 
0.1 < A < 1 and ^/i/vA in the observationally allowed range (110) or (111), both N^^^ and N^^^ vary only 
very little with the parameters (cf. figure 6). For the precision required here, it is enough to know N* at 
the precision of a whole number. Therefore, inserting parameters of the allowed order of magnitude, we can 
set the approximate values to 



min ^ ' 5 

N*„.a. = 59. 



(122) 
(123) 



In the two limiting cases, the spectral parameters can be evaluated at these values for N* . Neglecting the 
small variation of N* with the parameters of the theory, we observe that the amplitude P^ of the scalar 
spectrum depends on the combination ^/t/vA and on ^^. The expression(115) allows to understand the 
shape of the allowed parameter region in figure 3. The other spectral parameters, unlike the amplitude, are 
practically independent of ^^ and A and depend on the single parameter ^-^. Hence, we can plot n^, ac^ and 
r as functions of S,x (cf. figures 4 and 5). Compared to the numerical results, the approximate formula for 
the spectral tilt has an accuracy of the order of 10^'^ while the accuracy of the approximate results for a^ 
and r is of the order of 5 • 10^^. Given the uncertainties in the observational values, these accuracies are 
largely sufficient. 

As long as the quantity 4^^iV* is smaller than one (i.e. ^^ < 0.004), the series expansions of the hyperbolic 
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functions in (115)-(118) converge rapidly and the expressions can be further approximated by 

\/V*2 / -1 \ 

^c(fco) ^ y^^ (1 + 3(4ex^*)' + ■■■) , (124) 

Us (fco) - 1 ^ -il: fl + ^(4^x^*)' + •••) , (125) 



AT* V 3 
12 /_ 1 



"c(Ao) ^ -—2 ( 1 - ^i^^xN*f + ... ) , (126) 



KM :^ ^\^l--(MxNr + -j ■ (127) 

From these expressions, we can see that in the hmit ^^ — > 0, the predictions of the Higgs-Dilaton inflation 
model reduce to those found for the Higgs-Inflation model [31]. Hence, one can think of ^^ as the deviation 
of our predictions from those of the Higgs-Inflation model. In the Higgs-Dilaton scenario, the results for n^ , 
a^ and r in the limit ^t^, — > constitute a prediction of bounds on these quantities. We find (cf. figures 4 
and 5) 

n,(fco) < 0.97^1-^, (128) 

ac(fco) > -0.0006- ^, (129) 

r(fco) < 0.0033 ~ ^ . (130) 

These bounds are non-trivial predictions of our model. Further, given that Us, a^ and r are functions of ^^ 
only, the bound ^^ < 0.008 deduced from the observational lower bound on Ug translates to (cf. figure 5) 

ac(fco) < -0.00015, (131) 

riko) > 0.0009. (132) 

The upper bound on Us is well in accord with the observational constraints. Results of the Planck mission 
are expected to reduce the errors by a factor of a few and will hence provide an important test of the Higgs- 
Dilaton model [74]. While the present observational limits on a^ and r are too weak to compete with the 
bounds derived above ]64], the results of Planck might also improve this situation. 

We will see in section IV A that, if the scalar fields constitute a dark-energy component at late times, 
its equation of state parameter w^jg is also a function of the parameter ^^ only. As a consequence, the 
observational lower bound on Ug will induce a bound on w^j-,^. Thereby, the Higgs-Dilaton model provides 
a non-trivial connection between observables related to the early and the late universe. 

Besides the parameter region in which both ^^ and 4,^^N* are small, the observational bounds do not 
completely exclude the region where ^^ is small, but 4^^7V* is somewhat bigger than unity. In this region, 
as can be deduced from the approximate results (116)-(118), the predicted values for a^ and r go to zero 
exponentially with growing 4^^ TV*, while the spectral tilt becomes asymptotically linear in ^^ (cf. figure 4), 
i.e. 

n,(fco)-l^-8ex, for ^ < ^^ « 1 . (133) 

This fact will allow us to speculate about a somewhat deeper connection between n,, and wde in section 
IV A. 

From the approximate results (116)-(118) one can see that a^ and r are suppressed with respect to Us — 1. 
In the region where 4^^ TV* < 1, the suppression factor is 1/N*, while in the region where 4^^Af* > 1, the 
suppression factor is smaller than ^y^. This fact justifies the use of the observational bounds (104) and (105), 
that are based on the assumption of a power-law spectrum for scalar perturbations and the absence of tensor 
modes. 

The above results provide limits on the reheating temperature T^h, defined as the initial temperature of 
the homogeneous radiation dominated universe. Trh is related to Qrh through 

Qrh = 7^9cff{Trh)TX , (134) 
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where ges{Trh) is the effective number of relativistic degrees of freedom present in the thermal bath at 
the temperature Trh- Counting all degrees of freedom of the Standard Model plus the dilaton, one has 
ges{Trh) = 107.75. To lowest non trivial order in ^^ and 1/S,h, the limits on grh, i.e. ^(^X") < (^rh < V{6end) 
(cf. (108) and (109)) are found to be 

V{9Tn ^ T^Ml,, (135) 



A 



V{9,r,d) ^ -^XMl,, (136) 



where X = 7 — 4\/3 ~ 0.7. The bounds on Qrh, together with the obtained bounds on £,h/V^ (HO) and 
(111), translate into the following bounds on the reheating temperature T^h 




3.8 J^^ ■ 10^2 GcV < TrH < 2.5 • 10^'^ GeV . (137) 

V u 

Finally, let us note that the findings of this section allow us to constrain the region of initial conditions 
for the scalar fields which lead to successful inflation. Based on the assumption that the last N* e-folds 
of inflation take place in the scale-invariant region, we have found the fleld value 9* close to which the 
observable scales exit the Hubble horizon during inflation. The initial conditions for inflation have to be 
such that 9in > 9* . In terms of the original variables, this condition reads 



tanr . (138) 



6* was found to be 9* ~ arccos (e"'*^'^^ ) (114). For typical parameter values ^^ = 0.005, ^h — 65000 and 
N* — 58 one obtains 9* ~ 1.25 and -^^^^^ > 0.005. Considerations related to dark energy (section IV B) will 
yield an additional constraint on the initial conditions. The region of acceptable initial conditions satisfying 
both constraints is shown in Fig. 9. 

IV. IMPLICATIONS FOR THE LATE UNIVERSE 

We have shown in the previous section that a number of parameters of the theory can be constrained by 
two independent inflationary observables: the amplitude and the tilt of the primordial spectrum of scalar 
perturbations Vq, cf. Figure 3. The theory is therefore completely specifled at the inflationary stage, and 
any subsequent period should be consistent with that choice of parameters. In this section we focus on the 
late dark energy dominated stage previously described in section TIC, during which the dilaton fleld is rolling 
down along one of the potential valleys. In subsection IV A we show how this results in the dilaton playing 
the role of a quintessence fleld. We then derive consistency conditions among the inflationary observables 
and those associated to the dark energy dominated stage, which could allow us to either conflrm or exclude 
the model in the coming years. In subsection IV B, we derive a constraint that has to be imposed on the 
initial conditions of the scalar fields in order to have a successful description of dark energy. This will prove 
a posteriori that the whole period of observable inflation must take place in the scale-invariant region. 

A. The dilaton as quintessence field 

After the phase of reheating, the system enters the radiation dominated stage, at the beginning of which the 
total energy density is given by Qrh (cf. (134)). At that moment, the scalar flelds have almost settled down 
in one of the potential valleys, i.e. hit)"^ ~ f x(i)^ or, in terms of the variables {p,9), tanh {a9{t)/Mp) ~ 
1 — <; + 0{a). We will work with the assumption that the equality is exact and that the flelds 



1--; 



1+f 



i+etfe 



A i+esx 



28 

evolve exactly along the valley (cf. [1])^*'. In this case, at the level of homogeneous fields, we are left with 
a single degree of freedom p{t). As discussed in section II, also at the level of perturbations the field p 
(or equivalently p) is almost decoupled from the SM fields. Hence, we will from now on treat p as a field 
minimally coupled to gravity and not interacting with matter and radiation. Given that /i <C 1 and a ^ 1, 
its dynamics is described by the E-frame Lagrangian (25) with K and U = V + Vaq given by (50), (51) and 
(52), where one inserts the constraint tanh^(a^(i)/M) ~ 1 — i^, i.e. 

r /\/f2 ^ -1 

-j^ ^ -^R 7,{dpf - Vqe{p) , (139) 

with 

yQE{p) = ^e-'^'^^''- . (140) 

7 

As already pointed out in [1], it is remarkable that the exponential potential, which was proposed for 
Quintessence (QE) long time ago [34, 35, 75], appears automatically in the present model. It hence turns 
out that the dilaton field p can play the role of QE. 

Let us now discuss in more detail the infiuence of the field p on standard homogeneous cosmology. The 
equation of motion for the homogeneous field p = p{t) in spatially flat FLRW spacetime is given by 

p + 3Hp+^^=0. (141) 

dp 

Defining energy density qqe, pressure pqe and equation of state parameter wqe of the scalar field p as 

BQE = Ip^ + Vqe, (142) 

PQE = Ip^-Vqe, (143) 

WQE ^ ^^, (144) 

Qqe 

its equation of motion (141) can equivalently be written as 

QQE = -iHgQE{l + WQE) ■ (145) 

On the other hand, in the presence of a barotropic fiuid of energy density gi,, for instance relativistic or 
non-relativistic matter, the Hubble parameter is given by the first Friedmann equation as 

^' - ^ U + qqe) , (146) 

which in terms of the relative abundances fl = g/3MpH^ can be written as the cosmic sum rule 
ilb + ^QE = 1- The cosmological model described by equations (141) and (146) with a scalar field evolving 
in an exponential potential has been widely studied in the literature (for a recent review see [76]). We want 
to recap the main results established in the literature and show how they apply to our model. 

For the qualitative analysis of the system we rewrite equations (141) and (146) in terms of the observable 
quantities ^Iqe and Sqe = 1 + wqe as (cf. e.g. ]77]) 

S'qe = -3<5Qij(2 - 6qe) + 47(2 - SQE)^3SQEnQE , (147) 

n'QE = 3is,-SQE)nQEii-nQE), (i48) 

where prime, as before, denotes the derivative with respect to the number of e-folds N = In a. Further, 
i5b = 1 + Wb, where Wb is the equation of state parameter of the barotropic fluid. For radiation one has 



^^ The validity of this approximation can be checked numerically. Let us further note that the trajectory going exactly along 
the valley is an asymptotic, but not an exact solution of the equations of motion. 
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5b = 4/3, while for non-relativistic matter 5b = 1- An additional dark energy component with constant 
equation of state would have 5b < 2/3. In the scale-invariant model analyzed here, a component of this type 
is present as soon as the action (11) contains a term fix'^ with /3 > (cf. also section II A). For the reasons 
mentioned in section II A 3, we will mainly focus on the case /3 = 0, in which p will be alone responsible for 
dark energy.^^ It has been shown in [34, 36, 78] that as long a,s < 5b < 2 and depending on the value of 
the parameter 7, the system approaches one of two qualitatively very different attractor solutions. 

For 47 > ^/35b, the variables evolve towards the stable fixed point ^qe = 3(5^/167^ and 5qe — 5b- This 
means that the scalar field inherits the equation of state parameter of the barotropic fluid. Hence, the 
energy density of the scalar field scales like the energy density of the fluid. Unless the scalar field gives the 
dominating contribution to the energy density from the very beginning, it will never become dominating. 
Therefore, these so-called "scaling solutions" can not be responsible for the late-time acceleration of the 
universe. In this case, the accelerated expansion must be due to another mechanism, e.g. a barotropic dark 
energy component with 5b < 2/3. In other words, a scaling field can at best provide a small contribution to 
dark energy. 

For 47 < y/35b, the situation is very different. The stable fixed point is given by ^qe — 1 ^^nd 
5qe = 167^/3. Hence, in this case, the asymptotic solution describes a scalar field dominated universe, 
which is accelerating if 47 < v2, i-e. Cx ^ 5- This means that the scalar field with exponential potential 
and 47 < \/35b can describe the late-time acceleration of the universe, provided that the system has not 
quite reached the fixed-point by today^^. 

In the previous section, we have found that our model can successfully describe inflation if Cx ^ ^ ^ 10^"^. 
This yields the bound 47 ~ 4->/^ < 0.36. Hence, for this parameter choice, he system evolves towards 
the second type of fixed point, corresponding to a scalar- field dominated universe in accelerated expansion. 
This allows us to draw a non-trivial conclusion. Namely, if the parameters of the model are fixed by the 
requirements of inflation, and for /3 = 0, the late time behavior of the system necessarily corresponds to an 
accelerating universe, dominated by p. 

Current observations [64] show that the present abundance of dark energy is fi^j^; — 0.74. For /3 = 0, dark 
energy is entirely due to p and we can identify ^1%^ — f^og and Wq^ = w'jj^. The observed value shows 
that dark energy is not clearly dominating the present universe, which means that the system must not have 
reached its fixed point yet. 

We now want to qualitatively discuss the scenario in which the field p is irrelevant during the radiation 
and matter dominated stages, but has become important recently and is now responsible for the present 
acceleration of the universe. During the radiation and matter dominated stages one must have ^Iqe ^ 1- 
As long as this is the case, the second term on the right-hand side of (147) is small compared to the first 
one. Hence, 5qe is driven towards a very small value 5qe <C 1 and wqe — —1. This shows that even 
if initially qqe were dominated by kinetic energy, the kinetic part would soon die away and ^qe become 
potential dominated ^^ . As a consequence, the value of p is almost constant during the radiation and matter 
dominated epochs and remains practically equal to its value at the end of reheating. However, since qqe 
decreases more slowly than the energy densities of radiation and matter, ^Iqe becomes relevant at some 
point. At this point, p starts rolling down the potential and 5qe starts growing towards its attractor value. 
The initial value of ^qe has to be small enough such that Qqe remains negligible throughout radiation 
domination and only becomes important in the late matter dominated stage. The described scenario in 
which the quintessence field remains constant for a long time and then starts rolling down the potential 
goes under the name of "thawing quintessence" [37] . Two recent studies treating the case of an exponential 
potential can be found in ]77] and ]80]. 

In the approximation where 5qe = I+w^qb ^ 1, the system of equations (147) and (148) can be integrated, 
and one finds the interesting relation ]77] 

167^ 
l + WQE=^^F{nQE), (149) 



^^ Let us mention, that even in the case /3 < 0, appropriate initial conditions lead to a dark energy dominated phase. 



As has been shown in [79] , this statement holds even if 47 > v^- 
In principle one could imagine a scenario in which after reheating ^qe is non-negligible, as long as Sqe — 1- Since the 
kinetic part of Qqe decreases as a~^ it would soon fall below (Jjadiation and radiation would start dominating provided that 
the potential energy of p is small enough. However, we do not expect this to happen in our model, because the field p is 
almost constant during reheating. 
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Figure 7: The function F{Qqe)- Note that it becomes exactly 1/2 for JIqe = 0.739, which is very close to the 
observed abundance of dark energy Q.%e = 0.725 ± 0.048 ( WMAP7+BAO+Ho at 99% confidence level [64]). At this 
value the functional relation between the spectral tilt Us of CMB anisotropics and the equation of state parameter 
Wqe is particularly simple. 
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The function F{Q,qe) is monotonically increasing from to 1 (cf. figure 7). For the observed value ^q^ 



^''DE 



0.74, one gets F{n^Q^ == 0.74) ~ 0.5. Inserting this into (149) and identifying Wq^ 



-"DE' 



we obtain 



the following result for the present equation of state parameter of dark energy (provided by quintessence 
field p) 



1 



'-"DE 



8 Cx 
3 1 + ecx 



(151) 



We can now plug into this relation the upper bound ^^ < 0.008 (cf. eq. (Ill)), derived from the observational 
bound on Ug, as well as the theoretical lower bound ^^ > 0, such that (cf. Fig. 8) 



0< 1 



-"DE 



< 0.02 



(152) 



Thus, we have found that the parameter bound from inflation implies a strong bound on the equation of 
state parameter of dark energy. This is a rather non-trivial result. The current observational constraint 
-0.52 < 1 + w%j^ < 0.32 (WMAP7 + BAO + Ho at 99% confidence level [64]) is much too weak to 
compete with this theoretical prediction. From this point of view, the energy density qqe is practically 
indistinguishable from a cosmological constant. Nevertheless, the observational bound is expected to improve 
considerably in the near future. While the Dark Energy Survey collaboration aims at measurement of w^^ 
with an accuracy of ~ 5% [81[, the expected accuracy from the Euclid consortium is ~ 2% [82[. These 
measurements, together with the projected improvement on the determination of n^ from the Planck mission 
[74[, should provide an important consistency check of the Higgs-Dilaton model in the near future. 

In fact, the theoretical predictions of our model (with /3 = 0) can be further refined. Namely, since both 
the scalar spectral index Us and the equation of state parameter w'^^ depend mainly on ^^, it is possible 
to establish a functional relation between these two very different observables. Combining (151) with the 
approximate relation (116) allows us to express the scalar tilt n^ as a function of S'jje and the number of 
e-folds N* as 



^^S%E 

4 - 9S%^ 



coth 






(153) 



We plot this relation and the corresponding numerical result in figure 8. The plot is equivalent to the plot of 
Fig. 4, except that the independent variable is changed from £^^ to 5de = 1 + w'^^ with the help of (151). 
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Figure 8: This plot shows the approximate functional relationship between Us and w%e- The plain curves are 
numerical results. The red plain curve is obtained for Qrh = drh^ (instantaneous reheating), while the blue plain 
curve represents the case q^h = f?™" (long reheating). The dashed curves are obtained from the approximate relation 
(153). The red one for A'^* = N^ax = 59 (instantaneous reheating) and the blue one for A^* = N^^^ = 57 (long 
reheating). The black dashed straight line corresponds to the relation (155), approximately valid in the region 
j^ < 1 + w^iE ^ 1- Note that it does not exactly correspond to the straight line in figure 4, which would not 
correspond to a straight line in the present plot. The shaded region represents the experimental mean value and the 
associated la and 2a confidence intervals. 



As before, we see that the result is rather insensitive to variations of N* in the range A^^jj^ < N* < N* 
cf. (122) and (123). One can also derive a relation involving the respective second order quantities 



a^ 



48wf,E 



(-4 + 9<50,^)3 



(_4 + 9<5°,s)coth' ^^^ ^^^ 



4 — Qr?" 



6iV*(5^£;sinh" 



4 - %l^ 



(154) 



connecting the running uq of the scalar spectral index to the equation of state parameter wj^^ of dark energy 
and its rate of change wf-,^, defined through 

WDsia) = w%E + Woe ln(a/ao) . 

While the first order consistency condition (153) should become testable in the near future, a the test of the 
second order relation (154) will be more challenging. 

We have mentioned previously that the parameter region where rig is well-approximated by the asymptotic 
linear in £^^ is not excluded by observations. In terms of 5de this region is given by 5de <C 1 and 1 < 

3 



6Af gg) 
4_qAU 



^N*6l)^, in which the relation (153) becomes approximately 

2 



3(1 + u;^s) ~ (n, - 1) 



for 



37V* 



< 5de < 1 



(155) 



which can equivalently be written as a relation between "first orders" in the early and the late universe 

d In g^^ d In Pq 



rflnf 



cflnfc 



(156) 



Whether this is a fundamental consequence of SI or just a coincidence remains yet unclear. Note that if 
relation (156) should hold, it would not only imply that the deviation S'^^ of dark energy from a cosmological 
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constant is proportional to the deviation Us of the primordial spectrum from the scale-invariant one. In fact, 
if we could take it at face value, it would imply a concrete value for the present abundance of dark energy: 
F{^de) = 1/2 =^ ^DE = 0.739, surprisingly close to the observed value. In the same region of parameter 
space, also the respective second order quantities are proportional to each other, again for Que — 0.74, 

'Sw^DE - "C . (157) 

or equivalently 



(dlna)2 ^ (rfln/c)2 



(158) 



Let us stress again that the links between the observables Ug and a^, related to inflation, and w'^^ and 
w'j^Ei related to dark energy, are non-trivial predictions of the present model. They relate two a priori 
totally independent periods and allow us to use the measurable observables from CMB anisotropics to make 
predictions for the widely unknown DE sector. On the other hand, one should bear in mind that these results 
rely on several important assumptions. In particular, the functional relations are based on the requirement 
that the J-frame potential has a flat direction (/3 = 0) . 

B. Dark energy constraints on initial conditions 

Let us now show how the obtained results justify the assumption we made about inflation taking place in 
the scale-invariant region. From (152) we infer that qqe is dominated by the potential energy contribution 
and hence qqe — VQsip)- This allows us to deduce from the observational value Qq^ = ^ds — ^-'^'^ today's 
value of j5, 

1 7i<^ i__ / 4Aoff 



Po:.--Mpln^^^^j, (159) 

where we have defined an effective cosmological constant as 

Aoff = SM^H^n^E - 10"^^°M|, . (160) 

Now, numerical simulations show that the field p has been almost constant from the end of infiation till today. 
Therefore, the value of po provides an order of magnitude estimate for the value of p at the end of inflation. 
During the analysis of inflation we have made the assumption that the whole period of observable inflation, 
i.e. the last ~ 60 e-folds, took place in the scale-invariant region, where fi,f2 «C 1 and hence p ~ est.. We 
can now check this assumption by computing vi and V2, cf. (53) and (54), at p ~ p* ~ Pend — IPo- Using 
(159) and working in the usual approximation ^^ <C 1 and ^/^ 3> 1, we obtain 

144e|fAeff^ ,,24C,C^Aeff 1 

'" A M^sin^e' ^^- A Mlsin^^cos^' ^ ' 

From (112) and (114) we have ^end — 2 * 3i \/^ and 9* ~ arccos (e^^^>f^ ) . Evaluating vi and V2 for values 
^^, ^h and A^* of the orders of magnitude found in section HID and 0.1 < A < 1, we find that for the whole 
interval Oend < < 6*, ui,U2 ^tC 1, and hence that the deviation from exact scale invariance is negligible. 
This justifies a posteriori the neglecting of Ao during inflation. Let us note that this conclusion is not altered 
if one takes into account the slight change of the scalar flelds between the end of inflation and today. The 
change of p during the reheating oscillations and during the thawing quintessence stage are of the percent 
level. 

In section HID we have seen that for a successful description of inflation, the initial conditions for the 

scalar fields need to satisfy 6^ > 0* , respectively -^ > y i_|_gf tan6>*. We recall that for typical values 
Cx = 0.005, £,h = 65000 and N* = 58 one obtains ^^ > 0.005. The observational value for n%j^ ~ 0.74 
(respectively (159)) together with the knowledge that the field p remains almost constant from horizon 
crossing during infiation until today, allows us to further restrict the region of allowed initial conditions (cf. 
Fig. 9). Namely, if p is alone responsible for dark energy, and as long as the initial conditions lie in the 
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scale-invariant region (wi,'L'2 <C 1), the relation (159) yields approximately the initial value for the field p, 
i.e. 

Mp ( 4Aeff\ 

In terms of the original variables, this corresponds to a relation between Xin and hin given by (^^^ <C 1 and 
a » 1) 

-^+6a-i7^^7-TI7^-^10 . (163) 

Together with the bound /im/Xm ^ 10^'', this shows that initial conditions have to approximately satisfy 

hin/ ^Q ^ 10'^". Hence, the initial value of h has to be much larger than the scale Aq' . For p to exactly 
produce the observed abundance of dark energy, the initial values have to be chosen precisely on a line in 
the (p, 9)-, respectively the (x, ft.)-plane. This tuning of initial conditions is commonly referred to as the 
Cosmic Coincidence Problem (see e.g. [83])^^. Our model does not alleviate this problem with respect to 
other quintessence models. In fact, if one allows for an additional dark energy component, present if /3 > 0, 
the set of acceptable initial conditions extends to an infinite region. In that case, while the fine-tuning issue 
does not concern the initial conditions, the parameter /3 has to be finely tuned. Hence, in either case some 
"fine-tuning" is needed. At this point, it should be recalled that although the Cosmic Coincidence Problem 
is an undesirable feature, it is not a consistency problem and therefore does not invalidate this and other 
models of dynamical dark energy. 

Finally, we can briefly comment on the case of initial conditions lying in the region where Aq can not be 
neglected. Initial conditions lie in this region (w2 > 1) whenever 9in is sufficiently close to tt/2, respectively 
when hin/xin is sufficiently big. Note, however, that as a consequence of condition (162) and for typical 
parameter values, this only happens for extreme values 7r/2 — 6in ^ 10^^^^, respectively hin/xin ^ 10^°^. In 
that region p is no longer constant. The E- frame potential (29) becomes dominated by the term proportional 

to Aq, i.e. U{h,x) — yAo{h,x) — tt — 2!^ °h2\2 ■ The effect of this potential is to drive the scalar fields to 
larger values of x and h, respectively larger values of p, before they enter into the scale-invariant region. 
Qualitatively, this means that in the non-scale-invariant region the line of successful initial conditions is no 
longer given by (162) (respectively (163)) but turns towards the origin. Still, the discussion related to the 
Cosmic Coincidence Problem equally applies to initial conditions in this region. 

V. CONCLUSIONS 

We have considered a minimal scale-invariant extension of the Standard Model, non-minimally coupled to 
gravity, including a scalar dilaton. All mass scales at the classical level, including the Planck scale and the 
Electroweak scale, are induced by the spontaneous breaking of the scale invariance. The physical dilaton 
is almost massless but hardly affects particle physics phenomenology. Our findings rely on SI both at the 
classical and the quantum level [2, 3]. The replacement of standard General Relativity by Unimodular 
Gravity gives rise to an arbitrary constant in the equations of motion, which in a minimally coupled theory 
would play the role of a cosmological constant. However, due to the non-minimal couplings between the 
scalar and the gravitational sectors, this constant gives rise to a non-trivial "run-away" potential for the 
dilaton. As a consequence, the dilaton can play the role of a quintessence field, responsible for a late dark 
energy dominated stage. For appropriate values of the free parameters and initial conditions, the constructed 
model presents a rich cosmological phenomenology, providing mechanisms both for inflation and dark energy. 

We find that the amplitude P,^ of CMB anisotropics depends mainly on the ratio £,h/V^ while the spectral 
tilt Us, the associated running a^ as well as the scalar-to-tensor ratio r depend mainly on ^^. The observa- 
tional limits on Pq and n^ put bounds on ^/i/vA and S,x, which in turn provide the bounds a^ < —0.00015 
and r> 0.0009. In addition, the model predicts the bounds Ug < 0.97, a^ > —0.0006 and r < 0.0033, which 



For a discussion of the fine-tuning issue in the particular case of a Quintessence field with an exponential potential, see [79]. 
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Figure 9: This plot shows the different regions of initial conditions giving rise to qualitatively different evolutions. 
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For a successful description of inflation initial conditions have to lie above the line 9 — 0* c^ arccos ( e 

Ao > 0, the scalar fields contribute to dark energy in the late stage. Initial conditions have to lie above the arc 

of an ellipse given by p ~ 7po — -f- In (7*^^'^ ), for this contribution not to exceed the observed value of JIde- 

Hence, the yellow region corresponds to initial conditions giving rise to successful inflation and a contribution to dark 
energy not exceeding 0.jy^ = 0.74. The red segment of the ellipse corresponds to initial conditions for which the 
scalar fields yield the total observed dark energy. The hyperbola is given by f 2 = 1. Initial conditions below the 
hyperbola lie in the non-scale-invariant region, where Ao is important. Trajectories starting here tend to move away 
from the origin before entering the scale-invariant region and following a scale-invariant trajectory. Therefore, such 
initial conditions can also be acceptable as long as the corresponding trajectories enter the scale-invariant region at 
or above the line given by p. Note that while we only describe the quadrant x/^o ' ^/^o > 0, the reasoning would 
be completely analog in the other quadrants. 



are obtained in the limit ^^ — > and correspond to the predictions of the Higgs-Inflation model of [1]. The 
confrontation of these bounds with the results of the Plack satellite mission will constitute an important test 
of the Higgs-Dilaton model. 

Neither SI nor Unimodular Gravity forbid the existence of a quartic term Px^ in the Jordan frame, 
which would correspond to a proper cosmological constant in the Einstein frame. However, the parameter 
choice forbidding such a term (/3 = 0) appears to be specially interesting, both from the cosmological and 
the quantum theory point of view. For this choice, the dilaton is alone responsible for dark energy. The 
associated equation of state parameter w^,^ is found to practically depend on ^^ only. This has the interesting 



consequence that the spectral index n^ can be expressed as ris 



^DE. 



thus relating an observable from 



the very early universe to an observable of the present universe. For a particular parameter region, this 
relation takes the simple form — 3(w^^ -f 1) ~ {ug — 1). The observational bound on Ug translates into a 
bound < 1 + w'jj^ < 0.02, which might be tested by future experiments. Further, we were able to derive a 
relation between the running of the spectral index a^ and the rate of change Sw'j^^ of the equation of state 
parameter. Notice that for the dilaton to provide the measured abundance of dark energy, initial conditions 
have to be finely tuned. Hence, as is the case for all quintessence models, the Cosmological Coincidence 
Problem remains unsolved. 
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Appendix A: Higgs-Dilaton Inflation in the Jordan Frame 

This appendix is devoted to the study of the inflationary trajectories in the Jordan frame. We perform 
an analytical study of the trajectories of the scalar fields during slow-roll and compare it with the results 
of an exact numerical computation in Jordan and Einstein representations. The numerical computation in 
the Einstein frame is performed in the way described in [55] and takes into account the non-minimal kinetic 
terms in (25). As the classical level different frames just correspond to different choice of variables, and 
therefore the final physical results should not differ. Notice however the different units used in every frame. 
For homogeneous fields, h = h{t) and x — x(i)i the conformal transformation of the metric from the Jordan 
to the Einstein frame (22) depends only on time il = il{t) and implies a redefinition of the cosmic time, 
dt — il{t)dt , as well as the scale factor, a(i) = ^{t)a{t) , in the Einstein frame. The associated Hubble rate 
should be also redefined as 

where the prime denotes derivatives wrt the number of e-folds N in the Jordan frame. Relation (Al) allows 
us to easily obtain a useful relation between the number of e-folds computed in both frames, 

dN _, d\nn ,,„, 

A = — ^ = 1 ^ . (A2) 

dN dN 

Integrating this equation from the initial field configuration (f>o at the beginning of inflation we get 

N-N^\n§^<0. (A3) 

As expected, the number of e-folds is not an invariant under conformal transformations. However, the 
difference between the two frames turns out to be practically irrelevant during the inflationary stage. To 
obtain an upper bound, we focus on the value at the end of inflation, \nfl^nd/^Oi where the discrepancy 
between N and N is larger. As we saw in Section HI A 2, the inflationary are well described by ellipses with 
constant radius, rg = (1 + 6^h)^o + (1 + ^^x)Xo- Here ho and xo are the initial values for the Higgs and 
dilaton flelds respectively. Let us assume that they are roughly equal. In this case, it is possible to relate 
the initial and flnal amplitude of the h fleld to obtain 

(A4) 

where we have used ^h 3> ^^ ^ ^^^^ ^® the approximate relation among the field amplitudes at the end of 
infiation x — \ s^- Taking into account (A4), we obtain 17end/f^o — \ i+i2f ' which corresponds, for a 
typical value ^^ = 0.005, to 




N ~N 

N 



< 2% . (A5) 



It can be shown numerically that during most of the inflationary stage the difference among the number 
of e-folds deflned in both frames (A2), is indeed quite smaller than the previous bound. Given the small 
difference between the number of e-folds deflned in Jordan and Einstein frames, we will from now on identify 
N — N. Let us now consider the Higgs-Dilaton Lagrangian density in the Jordan frame (19). We will 
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Figure 10: Evolution ol the angular variable 2 as a function of the number of e-folds A*' and detailed view of the last 
60 e-folds. The green (dot-dashed) hues represent the approximate slow-roll solutions given by (A15), while the red 
(solid) and blue (dashed) curves correspond to the result of an exact numerical computation performed in the Jordan 
and Einstein frames respectively. 




Figure 11: Evolution of the Higgs h and dilaton x fields as a function of the number of e-folds A''. The green (dot- 
dashed) lines represent the approximate slow-roll solutions given by (A16), while the red (solid) and blue (dashed) 
curves are exact numerical results in the Jordan and Einstein frames respectively. 



assume that the initial values of the fields are such that they evolve within the scale-invariant region, in 
which the Ao-terni in (21) can be neglected. Far away from the valleys of the potential the contribution of 
terms proportional to a ~ 0(10"'^'^) can also be safely ignored. The Klein-Gordon equations of motion for 
homogeneous scalar fields are then given by 



3i?(/)''4-V^'''-i/^"i? = 0, 



(A6) 



where the Ricci scalar R for a FLRW geometry and is given by i? = 6{H + 2H^). In the Jordan frame 
field-space indices are raised and lowered with the Euclidean metric 6ij. The two Friedmann equations can 
be written as 



3HV(0) = Irk + V{^) - 3iJ9o/(0) : 



(A7) 



/(0)i? = -3(92 + 3Hdo)fi(b) - rk + 4T/(0) . 



(A8) 



If we assume the fields to be homogenous during infiation. together with the standard slow-roll approximation, 
k'Pa ^ V, (j>°' -^ V'"" and (\)°- <C i/0° the equations of motion for the scalar fields (A6), expressed in terms 
of the number of e-folds iV, become 



37J^ 



-^'" + ^/'"i?, 



(A9) 
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while the Friedmann equations (A7) and (A8) sinipHfy respectively to (note that / = Hf) 

V ~ 3H^ (/ + /') , (AlO) 

fRc^W - 9H^f . (All) 

In the last equation, we have assumed extended slow-roll conditions, namely 1 + 6£,a4'°' ^ ^{4') and 1 + 
Q^a<j>°' <C Hf{(t>), which should be checked numerically a posteriori. Equations (AlO) and (All) imply that 
the Ricci scalar can be approximated as i? ~ 12H^ (1 + /V(4/)), which does not correspond to the usual 
approximation H ^ H^ . Although it can be checked numerically that the contribution of the extra term 
f /{4:f) is indeed very small, it must be explicitly maintained to preserve the conservation of the dilatational 
current in the slow-roll approximation. Indeed, combining equations (A9), (AlO) and (All) we obtain the 
field space constraint 

{l + 6^^)xx' + a + Q^h)hh'^0, (A12) 

which, as shown in Section III A 2, gives rise (in those cases in which the Ao term can be neglected)to 
inflationary trajectories that can be properly described as ellipses in field space. Therefore, it will be useful 
to apply the same strategy of Section III A 2 and rewrite the problem in terms of polar coordinates (r, z) , 
defined as 



r^^(l + 6e.)/.^ + (l + 6e,)x^ z^W|^^t||^'^, (A13) 

where z — tan^, cf. (43). The evolution equation for the previous variables can be computed making use of 
(A9), (AlO) and (All) to obtain 



0' T^~4g. ,^7- 1 + -,), (A14) 



where (^ depends on the couplings ^h,^x ^-'^^ ^^^ given by (47). The previous equations can be easily solved 
to obtain the evolution of the radial and angular coordinates with the number of e-folds 

where rg and zq stand for the initial values of the fields. The comparison between the slow-roll solution (A15) 
for the z variable and the exact solutions obtained numerically in Jordan and Einstein frames is shown in 
Fig. 10. Notice that, as pointed out above, we have identified the number of e-folds computed in Jordan 
with that computed Einstein frame, iV ~ A^, given the small difference between the two during the whole 
inflationary period. As expected, the evolution of the dimensionless quantity z does not depend on the 
chosen frame. Making use of (A15) it is also possible to compute the corresponding values of the original 
Higgs and dilaton fields, which in terms of the z variable can be written as 

h{N) = ^^g= (1 + z-'{N)y'/' , x{N) = -^ig= (1 + z'iN))-'^' . (A16) 

The comparison with the numerical solutions is shown in Fig. 11. 
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